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SUMMARY 


This  final  technical  report  describes  the  result  of  two  programs  concerned 
with  high  power  laser  propagation  through  the  atmosphere  that  are  particularly 
important  to  the  Navy  high  energy  laser  program.  The  first  study  deals  with 
an  experimental  laboratory  simulation  of  moving  stagnation  zones  a problem 
which  can  be  quite  serious  for  the  Navy  application.  The  second  study  is  the 
propagation  and  hon-linear  interaction  of  a high  intensity  laser  with  water 
droplets  or  fog,  again  in  a laboratory  simulation  experiment  under  realistic 
conditions.  The  experimental  results  are  given  in  terms  of  non-dimensional 
parameters  which  can  be  interpreted  in  terms  of  real  outside  atmospheric 
conditions.  In  addition,  the  results  are  modeled  by  simplified  analytic 
expressions  which  should  find  wide  use  in  systems  analysis  of  Navy  scenarios.  ^ 

The  stagnation  zone,  a region  of  the  propagation  path  with  no  net  motion 
between  the  laser  beam  and  the  air  can  occur  in  a number  of  Navy  scenarios. 

It  is  partic\ilarly  serious  becaxase  in  such  a region  the  absorption  of  laser 
beam  energy  causes  more  severe  heating  and  thus  stronger  thermal  b"*  ~ ■'ming 
than  would  result  if  there  were  relative  motion.  Under  Navy  sponsox-ship  in 
a program  preceding  this  one,  we  studied  the  problem  of  a fixed  stagnation 
zone,  because  of  its  sin^jlicity  and  experimental  tractability.  We  determined 
that  the  initial  distortion  increased  in  severity  but  that  the  distortion  (and 
thus  intensity  on  target)  reached  a quasi-steady  state  value.  An  analytic 
model  was  found  which  predicted  quite  accurately  the  experimental  results. 

In  more  realistic  situations,  the  stagnation  zone  is  not  stationary  but 
has  an  axial  motion,  due  normally  to  changes  in  the  laser  beam  slew  rate. 

The  motion  of  the  stagnation  zone  can  influence  the  thermal  distortion  and 
the  influence  of  stagnation  zone  motion  is  the  subject  of  the  present  contract. 
Details  of  the  experiment  and  analytic  modeling  are  contained  in  Section  A. 
Basically  what  we  found  was,  under  realistic  conditions,  the  thermal  distor- 
tion caused  by  a moving  stagnation  zone  followed  the  quasi-steady  state 
distortion.  Thus  the  transient  behavior  of  the  intensity  on  target  is  des- 
cribed by  the  steady  state  solutions  of  the  stationery  zones  which  is  modeled 
analytically.  Thus  the  potentially  complicated  problem  is  reduced  to  a more 
simple  form  which  csui  be  handled  analytically.  One  other  aspect  of  this  prob- 
lem remains  to  be  investigated  and  this  is  the  effect  of  a noncoplanar 
engagement  geometry  where  the  problem  is  no  longer  confined  to  a two  dimensional 
phase  but  includes  a second  component  of  slew  due  to  a difference  in  altitude 
between  the  target  and  the  laser  source. 


The  second  area  of  this  finetl  report  discusses  the  results  of  the  ex- 
perimental program  to  examine  the  interaction  of  high  intensity  laser  radia- 
tion with  fog  with  particxilar  emphasis  on  the  non-lineEU"  thermal  blooming 
of  the  laser  beam  caused  by  the  water  droplet  absorption.  The  program  was 
prima.rily  experimental,  utilizing  an  8 kW  CW  CO^  laser,  a lU  meter  long  cell, 
seeded  with  micron  size  water  drops,  and  associated  diagnostic  equipment* 
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At  intensities  of  a kW/cm  and  scaled  Navy  scenario  conditions,  it  was  found 
that  the  lewer  interaction  led  to  vaporization  of  the  droplets.  This  vapor- 
ization had  two  positive  effects,  l)  it  reduced  the  linear  loss  of  laser  beam 
energy  due  to  partial  elimination  of  the  fog  and  2)  the  thermal  distortion  is 
significantly  less  than  that  predicted  by  the  unperturbed  conditions.  (Previous 
to  this  study  thermal  blooming  predictions  were  bsised  only  on  initieil  condi- 
tions). The  contribution  of  hydrosol  absorption  to  thermal  blooming  was  modeled 
by  an  analytic  expression  where  the  vaporization  of  particules  acted  as  a 
saturable  absorber.  Insufficient  data  was  obtained  under  the  present  funding 
to  verify  sJ.!  aspects  of  the  hydrosol  Induced  blooming.  The  details  of 
the  experimental  procedxire  and  results  are  contained  in  Section  B of  this 
report . 
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CW  THERMAL  BLOOJCNG  IN  A MOVING  STAGNATION  ZONE 

SECTION  AI 
INTRODUCTION 


This  report  describes  an  experimental  investigation  of  thermal  blooming 
occuring  in  stagnation  zones,  volumes  of  medium  in  which  there  occurs  no 
net  flow  transverse  to  the  laser  axis,  which  Taove  axially  between  the  tar- 
get and  the  laser.  The  study  focusses  first  on  the  expansion  of  the  data 
base  for  laser  beam  intensity  degradation  by  stationary  stagnation  zones 
over  a range  of  beam  focussing  ratios  and  attenuation  and  then  uses  this 
data  base  for  comparison  with  data  taken  with  moving  stagnation  zones. 

This  report  is  divided  into  five  sections:  Section  II  describes  the  apparatus 
and  some  of  the  experimental  results.  Section  III  reviews  the  theory  used 
for  understeuiding  the  effect  of  a stagnation  zone  on  laser  beam  propagation 
and  proposes  a modification  on  the  earlier  form  of  this  theory.  Section  IV 
presents  the  experimental  results  and  Section  V the  conclusions,  principeLLly 
that  thermal  blooming  in  a stagnation  zone,  stationary  or  moving,  is  charac- 
terized by  a single  thermal  blooming  parameter  with  the  temporal  variation 
of  the  cn-target  intensity  in  the  moving  zone  ceise  being  determined  by  the 
tine  dependence  of  this  parameter. 
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EXPERIMENTAL  APPARATUS  AND  INSTRUMENTION 


A.  INTRODUCTION 

The  apparatus  used  in  the  experiments  was  described  on3y  in  general  terms 
in  Ref.  1 and  several  of  its  properties  are  treated  in  more  detail  in  this 
section. 

B,  APPARATUS  KINEMATICS 

Tl.e  apparatus  is  shown  schematically  in  Fig.  Al.  The  absorption  cell  is 
mounted  on  e.  plate  which  is  driven  by  an  electric  motor-gearbox  aspeMbly  at 
speed  X . As  shown,  the  top  of  the  cell  is  free  to  rotate  about  tue  pivot 
point  af  xts  base  and  is  coupled  to  a second  motor-gearbox  assembly  which 
in  turn  is  attached  to  the  same  plate  on  which  the  cell  is  mounted.  This 
second  motor  drives  the  top  of  the  cell  at  a time  varying  speed,  x^^  + x^t, 
in  direction  opposite  to  x . The  way  in  which  x^  and  *3^  ere  controlled  is 
described  in  Ref.  1 as  is  ?he  way  in  which  the  speed  of  the  cell  top  is 
monitored.  For  the  situation  shown  in  '>^ig.  1,  the  position  y„  of  a stagna- 
tion zone  center  is  given  in  terms  of  the  cell  length,  I byCl) 


C.  LASER  OPTICAL  TRAIN 

The  optical  train  used  for  conditioning  the  laser  beam  in  th^se  experiments 
is  now  described  in  somewhat  more  detail  than  was  done  in  Ref.  1.  The  laser 
used  in  these  experiments  was  a water  cooled  electric  discharge  CO  laser 
approximately  two  meters  long  equipped  with  a 10  m radius  of  curvature  totally 
reflecting  end  mirror  euid  a 65  percent  reflecting  germanium  flat  decoupling 
mirror.  Zinc  selenide  or  sodium  chloride  flats  at  Brewster’s  angle  were  used 
intracavity  to  vacuum  seal  the  laser  discharge.  To  control  the  oscillation 
mode  of  the  laser,  an  adjustable  iris  wm  positioned  at  the  output  window 
and  the  iris  diameter  was  decreased  to  suppress  oscillation  at  modes  other  than 
TEM  . Care  had  to  be  taken  to  avoid  excessive  clcsure  of  this  Iris  as  it  was. 
founS  that  doing  so  apercured  the  gausslan  shaped  output  Intensity  spatiad 
distribution  of  the  laser  and  resulted  in  Fresnel  diffraction  and  spreading 
of  the  output.  The  largest  output  power  used  in  the  experiments  performed 
was  30  watts.  Losses  in  the  optical  train  were  found  to  increase  sli^tly  with 
time  BO  that,  for  30  watts  laser  output,  18  watts  were  typically  incident  on 
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the  absorption  cell  in  experiments  performed  early  in  the  study,  with  this  power 
decreasing  to  l6  watts  in  the  later  experiments.  An  adjustable  teDescope,  con- 
sisting of  two  germaniun  lenses  of  +6,8  cm  and  +13  cm  focal  length,  was  used  to 
focus  the  laser  output  onto  the  detector.  These  lenses  are  shown  in  Fig.  A1 
to  be  in  the  "beam  expander"  sequence  with  the  shorter  focal  length  lens 
positioned  first  in  the  beam.  With  this  configuration  u value  of  about  1»,5 
could  be  attained  for  the  degree  of  focus,  F,  defined  as  the  ratio  of  beam 
diameter  at  the  cell  entrance  to  the  diameter  at  the  cell  exit.  With  the 
lenses  interchanged,  the  value  of  F attained  was  1.2.  The  laser  was  turned 
on  and  off  in  these  experiments  by  withdrawing  and  inserting  a knif'>  edge 
shutter  at  the  beam  waist  in  the  telescope,  the  point  were  the  beam  diameter 
is  the  smallest  in  the  entire  optical  train.  The  time  elapsed  in  withdrawing 
the  knife  edge  from  the  beam  waist  was  determined  by  monitoring  the  laser  power 
at  the  detector,  with  the  pivoting  mirror  stationary,  as  a function  of  time 
while  the  knife  edge  was  withdrawn.  The  turn-on  time,  defined  as  the  time 
elapsed  for  the  detected  signal  to  reach  a fraction  equal  to  (l-e”^)  of  its 
fxill  value,  was  measured  to  be  0.6  msec  when  the  telescope  was  in  the  beam 
expander  mode.  This  "time"  is  of  course  determined  by  the  speed  of  the  knife 
edge  and  the  diameter  of  the  beam  at  the  waist,  the  latter  being  equal,  to  the 
product  of  the  focal  length  of  the  first  telescope  lens  and  the  laser  beam 
divergence  angle.  Accordingly,  the  turn  on  time  in  the  case  where  the  longer 
focal  length  lens  is  the  first  telescope  lens  is  seen  to  be  1.6  msec.  It  is 
important  that  the  turn  on  time  be  small  compared  the  time  characteristic 
for  CW  thermal  blooming,  defined  as  the  time  elapsed  for  the  peak  intensity 
at  the  detector  to  equal  e of  its  unperturbed  value.  For  the  most  extreme 
experimental  conditions  encountered  (l6  W laser  power,  ten  atmospheres  COg,  F » 
U,8)  the  blooming  time  is  calculated,  from  a formula  given  in  the  appendix,  for 
the  apparatus  to  be  a 8.8  msec  and  therefore  large  compared  to  the  turn  on  time 
so  that  the  finite  magnitude  of  the  latter  will  not  affect  the  experimental 
results.  At  the  lower,  F * 1.2,  focussing  ratio,  the  blooming  time  is  75  msec, 
and  the  same  conclusion  holds. 

The  focussing  ratio,  F,  as  will  be  seen,  is  of  importance  to  the  parameters 
which  characterize  thermal  blooming  severity.  For  a beam  of  gaussian  spatial 
intensity  distribution  the  beam  radius  is  defined  as  the  displacement  from  beam 
center  at  which  the  intensity  is  e”  of  its  value  at  center,  which  in  turn  can 
be  expressed  in  tepis  of  the  beam  diameter,  d*  and  the  total  beam  power  P by 
the  ratio  P/(rAd»  ).  Alternately,  the  value  of  d*  is  the  diameter  within  which 
a fraction  equal  to  (l-e*"  ) of  P is  contained.  The  value  of  P is  calculated  by 
first  measuring  the  value  of  d*  at  the  exit  of  the  telescope  and  at  the  detector 
onto  which  the  telescope  focus  the  beam  and  using  the  relationsh’  ..  which  expresses 
the  dependence  of  d«  at  distances  z from  the  detector: 


d*(z)»(d*  (0)*+  Z*d*(L)*  L A-n-2 
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where  L,  the  distance  from  the  detector  to  the  telescope,  is  the  focal  length 
of  the  latter. 

The  quantity  d*(L)  is  obtained  by  installing  an  adjustable  iris  at  the  out- 
put lens  of  the  telescope  and  measuring  P{d),  the  power  at  the  detector  when 
the  iris  diametergis  d.  The  relevant  plot  of  this  experimental  data  is 
ln(l-P(d)/P)  vs  d , and  d*{L)  is  the  value  of  d for  which  the  ordinate  has  a 
value  of  minus  one.  One  such  plot  is  shown  in  Fig.  A2  from  which  it  is  seen 
that  the  value  of  d«(L)  is  1.1  cm. 

The  value  of  d*(0)  is  obtained  using  the  pivoting  mirror  shown  in  Fig.  A1 
in  conjunction  with  the  detector,  liquid  Nitrogen  cooled  gold  doped  germanium 
equipped  with  a 0.02  cm  pinhole  aperture,  which  samples  the  beam  intensity  as 
it  is  swept  across  the  latter  by  the  mirror.  The  mirror  is  mounted  in  a general 
scanning  G300PD  motor  driven  by  a Tektronix  AF501  Band  Pass  Filter  run  as  an 
oscillator.  The  beam  turned  by  the  mirror  follows  an  oscillatory  path  and  is 
sweeped  across  the  detector  aperture  at  a speed,  2n\)B  d,  determined  by  the  dis- 
placement of  the  aperture  from  the  mirror,  d,  the  oscTllator  frequency,  \J,  and 
the  maximum  angular  deflection  of  the  mirror,  6 . The  detector  output, 
terminated  by  a 10 kfl resistor  (resulting  RC  response  time  is  1 psec)  is  ampli- 
fied by  a second  AF^Ol  Band  Pass  Filter  and  displayed  on  an  oscilloscope.  One 
such  oscillogram  is  presented  in  Fig.  A3.  The  value  of  d*(0)  is  given  by  the 
product  of  the  sweep  speed,  8|0  cm/sec,  and  the  time  intearval  for  which  the 
measured  intensity  exceeds  e”'^  of  its  maximum  value.  For  the  sweep  shown  in 
Pig.  A3,  the  value  of  d*(0)  is  seen  to  be  0.095  cm,  five  times  larger  than  the 
pinhole  and  which,  when  combined  with  d*(L)  the  positions  of  the  enirance  (i2l*  cm) 
auid  exit  (20  cm)  windows  in  Eq.  2 (A-II-2)  gives  a value  of  •♦.B  for  the  degree 
of  focus. 


D.  GAS  HAI^DLING 

One  of  the  quantities  to  which  thermal  blooming  is  sensitive  is  the  product 
of  the  linear  attenuation  coefficient,  «»  and  t,  the  length  of  the  range  or, 
in  laboratory  experiments,  the  gas  cell.  Therefore,  two  different  gas  mixtxires 
were  employed  in  the  present  experiments.  For  at  values  up  to  0.2  a mixture 
of  nitrogen  and  propylene  was  xxsed  while  for  a value  of  0.1*5*  pure  GO  was 
used.  For  either  mixture  the  total  gas  pressxare  was  ten  atmospheres  absolute 
to  minimize  the  role  of  conduction  relative  to.gonvection  in  determining  the 
theirmal  blooming  properties  in  the  experiment. 

E.  DATA  COLLECTION 

In  a stationary  medium,  the  effect  of  thermal  blooming  is  to  spread  the 
laser  power  symmetrically  about  the  beam  axis,  but  putting  the  siediuffl  in  motion 
causes  the  spreading  to  be  symmetric  only  about  the  plane  passing  through  the 
beam  axis  and  parallel  to  the  medixaa  velocity,  with  the  maximum  intensity  in 
this  plane  not  always  coinciding  with  the  xinperturbed  beam  axis.  Since  the  ^ 
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quantity  of  interest  in  the  present  studj'  was  the  behavior  of  this  maximum 
intensity,  the  data  records  consisted  of  sweeping  the  beam  across  the  detector 
«q)erture  in  the  plane  described  above  at  a frequency  of  l66  Hz  and  observing 
the  variation  of  the  maximum  intensity  in  this  plane  with  x^,  x^  and  x^^  for 
various  combinations  of  F,  at,  and  P. 

Figure  aU  presents  the  data  records  of  two  experimental  runs.  Each 
record  consists  of  two  traces.  The  top  trace  in  each  case  is  the  laser 
detector  output  as  a function  of  time.  This  output  appears  as  a series  of 
spikes,  each  of  which  represents  a sweep  of  the  beam  across  the  detector  face 
as  described  above.  The  bottom  trace  represents  the  output  of  a GE  H13B1 
Photon  Coupled  Interrupter  Module,  essentially  a compact  light  source  and 
detecoor,  which  measures  the  amount  of  light  at  a fixed  point  through  the 
photograph  of  a grid  of  alternating  zones  of  high  and  low  opacity.  This  grid, 
as  indicated  in  Fig.  Al,  is  fixed  to  a point  displaced  7^*  cm  from  the  cell 
pivot  so  that  x is  determined  by  noting  the  elapsed  time  between  two  maxima 
in  the  bottom  trace  and  combining  this  with  the  known  distance,  0.055  cm 
between  opacity  minima  on  the  film  and  t)  e knowledge  of  the  position  of  the 
film  relative  to  the  cell  pivot.  In  Fig.  AUa  x is  observed  to  be  zero  and 
Xj  has  a value  of  5.^  cm/sec.  Jn  contrast,  the  F"  g.  AUb  indicates  a value  of 
5.7  cm/sec  for  x^  and  27  cm/sec  for  x^.  For  both  experiments  x^  had  a value 
of  -li.O  cm/sec  so  that  Fig.  AUa  represents  thermal  blooming  in  a stationary 
stagnation  zone  situated  at  75  cm  above  the  entrance  window  and  in  Fig.  A^^b 
the  stagnation  zone,  situated  initially  69  cm  above  the  entrance,  moves  toward 
the  entrance  window  as  given  by  Eq.  A-lI-1.  It  is  seen  that  the  laser  inten- 
sity for  these  two  runs  shows  a different  temporal  development,  becoming 
constant  at  a relatively  low  value  in  Fig.  AUa  but  passing  through  a minimum 
and  increasing  markedly  in  Fig.  AJ»b.  The  reasons  for  this  temporal  develop- 
ment will  be  discussed  in  later  sections. 

Generally,  the  quantity  of  interest  for  comparison  with  analytic  models  of 
thermal  blooming  is  ot  the  measured  intensity  but  rather  the  intensity  norma- 
lized by  that  which  ould  be  incident  to  the  detector  in  the  absence  of  thermal 
blooming  or  the  intensity  as  indicated,  for  an  empty  cell,  by  the  maximum  in 
the  trace  in  Fia.  A3  corrected  only  for  the  linear  attenuation  due  to  molecular 
absox*ptlon,  e * . When  al  is  small  as  in  Fig.  U,  the  thermal  blooming  time  is 
large  enough  and  the  intensity  of  the  bloomed  beam  large  enough  to  permit  both 
the  temporally  varying  intensity,  l(t),  and  the  normalizing  intensity,  l(0)  to 
be  read  from  the  same  oscillogram.  For  larger  al  values  it  was  found  that  1(0) 
was  an  elusive  quantity  to  measure.  One  reason  for  this  was  that  the  thermal 
blooming  time,  on  the  order  of  8.8  msec  was  comparable  to  the  uncertainty 
in  the  zero  time  point  in  the  oscillogram  due  to  the  fact  that  the  oscillo- 
scope was  triggered  by  the  first  beam  sweep  across  the  detector  which  could 
lag  behind  the  true  beam  tum-on  at  the  knife  wedge  shutter  by  half  the  pivot 
mirror  sweep  period  or  6 msec.  In  addition,  for  large  al  the  intensity  of 
the  bloomed  beam  was  so  small  as  to  require  additional  aaqilification  at  the 
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oscilloscope  to  facilitate  its  measurement  in  which  case  l(.0)  was  off  scale 
on  the  oscillogram.  In  this  situation  empty  cell  inteasity  traces  such  as 
that  shown  in  Pig.  A3  were  taken  as  c^ibration  points  and  l(o)  was  then 
the  product  of  this  intensity  and  e~  and  could  be  rendered  useful  to 
normalizing  the  data  on  an  oscillogram  merely  by  multiplication  by  the  addi- 
tional amplification  factor  used  in  obtaining  the  oscillogram. 
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SECTION  AIII 

anaijYtical  models  for  thermal  blooming  by  stagnation  zones 


A.  INTRODUCTION 

As  past  studies  have  shown  , thermal  blooming  is  a phenomenon  which 
degrades  laser  propagation  and  is  influenced  by  a number  of  factors  which 
include  beam  power,  P,  the  degree  of  focusing,  P,  the  linear  attenuation 
over  the  range  ai,  and  medium  properties  such  as  density,  p,  heat  capac- 
ity, c , and  refractive  index  temperature  coefficient,  It  has  been 

found  Selpfxil  to  devise  models  for  thennal  blooming  under  the  varying  condi- 
tions of  interest  to  derive  thermal  blooming  parameters  which  include  in  a 
single  quantity  the  properties  listed  above  along  with  properties  peculiar  to 
the  condition  of  interest^  . For  the  stationary  stagnation  zone  problem, 
this  was  done  by  Berger  * who  modeled  the  stagnation  zone  as  a fixed  thin 
thermal  lens  of  constant  properties  which  distorted  the  on-target  intensity 
by  shifting  the  focal  point  of  the  beam.  In  this  section,  a modified  form  of 
Berger's  treatment  is  introduced  to  provide  an  alternate  parameter  for  charac- 
terizing thermal  blooming  in  the  presence  of  a stagnation  zone, 

Berger's  thin  thermal  lens,  treatment  strictly  holds  for  stagnation  zones 
only  for  those  cases  in  which  the  properties  and  position  of  the  lens  are 
constant  in  time.  In  an  attempt  to  obtain  a treatment  of  stagnation  zone 
thermal  blooming  which  includes  non-steady  state  effects,  a second  model,  also 
has  been  developed  and  is  described  in  Appendix  AI, 

B. *  STEADY  STATE  THIN  THERMAL  LENS  MODEL  OF  STAGNATION  ZONE  THERMAL  BLOOMING 

(U) 

As  stated  originadly  by  Berger  the  effect  of  a stagnation  zone  on  laser 
beam  propagation  can  be  seen  by  considering  a thin  thermal  lens  of  negative 
focal  length  to  be  positioned  in  the  beam  at  the  stagnation  zone.  The  focal 
length  of  this  lens  is  given  in  terms  of  the  medium  properties,  the  beam 
radius,  a(z'),  at  z',  the  location  of  the  stagnation  zon's  in  the  range,  the 
on-time  of  the  laser,  t,  and  the  thickness  of  the  lens,  6z,  a quantity  to  be 
dealt  with  further  at  a later  point; 


f=  ir^CpO(z')V(^^aPt  d?). 


(A-III-l) 


7 


^ ..i*  ,1..  ....I  ■«»>■■■.  . fiYinif  t !t 


K77-92257&-13 


The  unpertvirhed  laser  beam,  propagating  in  the  z direction,  is  considered  to 
have  a radixis  a at  its  origin,  to  converge  to  a focus  a distance  B’  beyond 
z* , and  the  effect  of  the  t'.  in  lens  is  to  increase  the  distance  to  a value  B 
whose  magnitude  is  given  by  the  thin  lens  lav: 

R'+z't\ 


t(Q(Az') 


where  R'f /(R'+f ). 

Of  interest  is  the  intensity  at  a point  Zj.  between  z'  and  B*+z*  where 


l(Zf)=P/irQ'(Zf)^), 


(A-IIl-2) 


and  a'(z)  denotes  the  z dependence  of  the  perturbed  beam  radius; 


q‘(z)*  OolR’/ReHRet  /(«'+  ?') . 


(A-111.3) 


The  variation  l(*j)  with  the  position  and  properties  of  the  lens  is  readily 
expressed  by  ' 


0gzf)  ^ -2 

Kzf) 


where  the  first  two  terns  in  brackets  relate  to  lens  position  while  the  third 
relates  to  lens  properties.  In  principle,  solution  to  the  moving  stagnation 
zone  problem  is  given  by  the  solution  of  Eq.  (A-III-I4).  In  practice  it  is 
convenient  first  to  fix  the  position  of  the  lens  and  consider  the  stationary 
stagnation  zone  problem.  In  this  case,  Eq.  (A-III-U)  is  readily  integrated  to 
give 
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I = 111l:!L  (\’  ^l£!!il“ 

^ l{z^,oJ  \ a(2fjT  j 'A-iii-5) 


At  this  point  in  his  development,  Berger  points  out  that  there  comes  & 
time  after  which  the  increase  in  engergj , Pt,  deposited  at  the  stagnation  7.one 
to  cause  the  lens*  focal  length  to  decrease  is  balanced  by  the  decrease  in 
ga  due  to  convective  cleanout  of  heated  medium  ly  motion  of  the  latter  normal 
to  the  beam  axis  at  points  other  than  the  r.tagnation  zon:  center.  Once  this 
steady  state  is  reached,  Berger  postulates  that  is  replaced  by  the  quantity 
Ua(z')/w,  where  w is  both  the  angtilar  velocity  of  the  medijm  about  the  stagna- 
tion zone  and  also  the  slew  rate  of  the  ttrget  about  the  laser.  Mating  use  of 

-az~ 

the  analysis  of  the  apparatus  in  Ref.  1 and  including  the  factor  e which 
Rccotints  for  the  attenuation  of  the  laser  beam  pwer  in  the  path  preceding  the 
stagnation  zone,  Eq.  A-IIT-5  can  be  written  in  terms  of  the  apparatiiis 
parameters : 


Irss*  (i  + Nq)-,^ 

Kzf.o) 


(A-IlI-6) 


where 


wpcp  a 


(A-IIT-7) 


In  practice,  ^ was  found  not  to  follow  strictly  Eq.  {A~III-6)  but  i^|-^ead 
to  be  character! zed  by  the  thermal  blooming  parameter  R according  to  ' 


Irss*  <•+  0.85  Nq+  0.09 Nq*)*'. 


(A-III-8) 


The  credibility  of  thlj  .wdel  for  stagnation  zone  thermal  blooming  was  e.nhanced 
further  by  the  finding  that,  in  addition  to  being  a useful  parameter  for 
characterizing  the  thermal  blooming  data,  N.  correctly  predicts  the  scaling  of 
thermal  blooming  severity  with  the  focus ing^atlo,  F,  this  scaling  being  linear 
in  F in  the  limit  of  large  F. 
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Central  to  Berger's  treatment  is  the  statement  that  t&s  equaling  Ua(z')/w 

expresses  the  steady  state  situation  that  convective  c],eanout  of  heated 

medium  at  the  ends  of  the  range  places  the  limit  to  the  strength  of  the  Ions; 

to  Eilternative  treatment  is  the  statement  that  tSz  simply  he  on  the  order  of 

R /w  or  a (z  )/w  after  the  establishment  of  the  steady  soate.  The  replacement. 

o?  tS.z  by  a vz^)/w  .restilts  in  a therijal  blooming  distortion  parameter  which. is 

Independent  of  F at  large  F,  in  variance  with  Fleck's  analytical  resvlts. 

In  contrast,  replacing  t5z  with  a /w  gives 

o 


„ ^ . (-/i  T ) O ^ 0 > CXq/X*  ) . 


(A-III-9) 


as  the  relevant  blooming  parameter  in  the  expression 


Irss  = (1+Nl)*^ 


(A-III-IO) 


It  is  fotod  that  Eq.  (A-III-IO),  like  Eq.  (A*ITI-6),  is  not  strictly  adhered 
to  but  instead  an  expression,  similar  to  {A-III-8)  is  developed  involving 
to  ercpress  the  stationary  stagnation  zone  data. 

At  this  point  is  is  interesting  to  consider  some  properties  of  Eq. 

( A-III-9).  Of  particxilar  interest  is  the  severity  of  thermal  blooming  and  .its 
dependence  on  stagnation  zone  position  and  degree  of  focus.  This  question  is 
simplified  by  considering  the  fvinction 


X(X/,)Co.F)  = 


(FX/+{F-|)Xo',^ 


(A-lII-11) 


which  contains  all  the  kinematic  and  focusing  dependence  of  . This  function 

is  plotted  vs  x,  in  Fig.  A5  for  the  indicated  values  of  x and^ol,  the  value  of 
& o 

the  latter  being  characteristic  for  a one  meter  cell  conta..ning  ten  atmospheres 

of  COg.  It  is  seen  that  the  effect  of  increasing  F is  to  caiose  x to  exhibit  a 

sharper  and  higher  maxim\i>n.  At  large  x^,  when  the  stagnation  zone  is  close  to 

the  laser,  the  value  of  is  seen  to  tend  to  be  linear  in  F.  ConverseJy,  for 

a stagnation  zone  positioned  progressively  closer  to  the  eource^  the  value  of 

I is  exp(ccted  to  exhibit  a minimum  with  a relatively  sharp  recovsiqr  for 

larger  values  of  F and  a shallower  minimum  for  smaller  F values.  As  seen  In 

Fig.  Al)b  the  vadue  of  ly  does  indeed  exhibit  a minlmim  followed  by  a marked 

recovery  as  the  stagnation  zone  moves  toward  th.  source,  suggesting  that  for 
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the  gSB  mixture  and  zone  velocity  considered,  the  instantaneous  value  of  Ip 

is  given  by  I ^ corresponding  to  the  ‘instantaneous  value  of  As  is  the 

case  for  iC^fs  not  expected  that  the  measvired  value  of  I to  be  given 
Q / ^\  rss 

in  terms  or  exactly  as  in  Eq.  (A-III~6)  but  that  I show  a functional 
Xi  i^ss 

dependence  on  the  many  quantities  of  importance  to  blooming  in  the  presence 

of  a stagnation  zone  through  the  single  parameter  N . 

L 
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SECTION  AIV 
EXPERIMENTAL  RESULTS 


A.  INTRODUCTION 

The  results  of  the  measurements  of  Ij.  for  propagation  through  a stagnation 
aone  are  presented  in  this  section.  In  the  first  subsection,  properties  of 
stationBOV  stagnation  zones  are  Investigated,  specifically  the  dependence  of 
If  on  linear  attenuation  over  the  range  and  degree  of  focus,  and  the  experi- 
mental results,  along  with  those  of  Berger(5)(o)  and  analytical  results  by 
several  workers  as  presented  by  Berger are  shown  as  functions  of  Nq  and 
N^.  In  the  second  subsection,  thermal  blooming  in  a moving  stagnation  zone 
is  examined.,  again  as  a function  of  attenuation,  degree  of  focus  and  laser 
power.  The  prinicple  conclusion  is  that  for  the  zone  velocities  considered, 
the  time  dependence  of  Ij.  is  determined  by  the  time  dependence  of  N^  or  Nq 
and  the  sensitivity  of  Ij.  in  the  stationary  zone  case  on  these  parameters. 

B.  STATIONARY  STAGNATION  ZONES 

As  was  discovered  by  Berger^^^  and  as  is  indicated  in  Pig.  AUa,  the  value 
of  Ij.  for  a beam  propagating  through  a stationary  zone  experiences  a transient 
decay  before  assuming  a constant  vsdue  whose  magnitude  is  determined  by  the 
properties  of  the  propagating  medium,  the  laser  beam  and  the  kinematics 
determining  the  stagnation  zone  position  in  the  range.  Data  obtained  from 
oscillograms  as  that  in  Fig.  Al+a  is  shown  in  Fig.  a6  which  illustrates  the 
dependence  of  Ij.  on  kinematics  for  the  indicated  values  of  focusing  degree 
and  attenuation  over  the  range.  Qvialitatively , the  value  of  Ij.  is  seen  to 
behave  in  a manner  consistent  with  that  shown  for  x(xjj^,  x^,  5)  in  Pig.  A5, 
exhibiting  a sharp  drop  as  the  net  wind  speed  at  the  exit  window  increases 
through  zero,  a minimum  soon  thereafter  and  gradual  increase  at  still  larger  net 
wind  speeds.  This  behavior  gives  rise  to  the  expectation  that  Njj  as  well  as 
Nq  is  a suitable  peurameter  characterizing  thermal  blooming  in  a stagnation 
zone. 


The  value  of  Ij.  is  plotted  as  a function  of  Nq  in  Pig.  A7.  Also  included 
is  experimental  data  taken  by  Berger ^5)  analytical  results  obtained  by 
Ulrich  and  Wallace  and  Fleck,  as  presented  by  Berger^^)^  ajid  Fleck'"^^.  The 
di.ta  taken  in  the  present  study  is  seen  to  be  self  consistant  for  a fixed 
value  of  P.  Por  Berger's  data  the  value  cf  Nq  was  calculated  from  Eq,.  16  in 
Ref.  5 assuming  a linear  taper  in  the  laser  beam  radius  with  range  position. 

In  addition  Berger's  function  R{x)  was  set  equal  to  unity  Independent  of  zone 
position. 
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In  comparing  Berger’s  data  with  the  present,  it  is  recalled  that  in 
some  of  his  experiments,  Berger  used  a detector  fixed  at  the  cell  exit  so  as 
to  sample  the  intensity  only  on  the  axis  of  the  unperturbed  beam^  while 
in  others  he  used  a rotating  mirror  to  sweep  the  beam  across  the  detector 
with  the  cell  positioned  hori2ontally  thereby  giving  :^se  to  bouyancy  generated 
gas  flow^^^  transverse  to  the  beam  axis.  Thus,  in  the  former  experiments, 
Berger  measured  a property,  the  "on-axis'*  intensity,  different  th*ui  that 
meaisured  in  the  latter  and  in  the  present  experiment,  the  "maximum  intensity". 
As  calc\ilated  by  Fleck 'T)  the  latter  is  expected  to  exceed  the  former  by  as 
much  as  a factor  of  two.  However,  it  is  seen  that  the  scatter  in  the  experi- 
mental results  obscures  this  difference  except  for  the  tendency  in  Bergers 
data  to  lie  slightly  but  consistently  lower  than  the  present  at  a given  value 
of  the  distortion  parameter. 

In  Fig.  a8  the  dependence  of  Ip  on  is  presented.  The  only  discemable 
difference  between  the  presentations  in  Figs.  AT  and  A8  is  the  tendency  of 
the  data  to  be  some.That  more  tightly  clustered  about  the  fi’_  ^ed  curve  when 
plotted  in  terms  of  Nj^.  Agx'eement  with  the  analytical  results  is  good  using 
either  Nq  or  Nl.  For  comparison,  Ij.  is  T«iotted  against  Nl' , obtained  using 
a(yf)  in  the  steady  state  expression  for  t^z,  in  Fig.  A9.  Presented  in 
terms  of  Nj^*,  tight  clustering  about  the  fitted  ciirve  is  again  observed  par- 
ticularly for  the  low  F value  data.  However,  agreement  with  Fleck's  calcula- 
tions for  large  F is  very  poor,  indicating  that  Nl'  exhibits  incorrect 
sensitivity  to  degree  of  focus. 

» 

The  stationary  stagnation  zone  data  presented  here  was  taken  to  expand 
the  data  base  used  in  studying  moving  stagnation  zones.  In  the  next  sub- 
section, moving  stagnation  zone  data  is  examined  and  compared  to  the  station- 
ary zone  data  presented  above. 

C.  MOVING  STAGNATION  ZONES 

As  seen  in  Fig.  AUb  the  value  of  Ip  for  the  moving  stagnation  zone  case 
may  exhibit  significantly  more  temporal  veuriation  than  does  its  value  in  the 
stationary  zone  case. 

Tlie  data  of  Fig.  A^b,  representing  one  of  several  data  runs  taken  at  the 

indicated  ai  and  F values  with  varying  values  of  y^  and  A..  (A.  = Xj^/x  ),  is 

shown  in  its  reduced  form  in  Fig.  AlO  along  with  information  obtained  »om 

steady  state  data.  The  data  point  denoted  by  a cross  is  the  measured  value 

taken  on  the  same  day,  of  the  steady  state  ly  in  a stationary  zone  experiment 

with  the  same  and  x^  as  occurred  at  the  instant  of  time  indicated  in  the 

Koving  zone  eaqperiment.  The  coincidence  between  this  point  and  the  moving 

tone  data  indicates  that  I^  for  a moving  zone  characterized  by  temporally 

varying  * and  x is  determined  at  any  instant  by  the  steady  state  value  of 
A o 
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ly  corresponding  to  the  values  of  x and  x at  that  instant.  The  solid  line 
indicates  the  value  of  calculated  from  ?he  cvirve  drawn  through  the  data 
in  Fig.  A8.  The  agreement  between  this  "vurve  and  the  data  is  particularly 
good  with  the  curve  tracing  the  temporal  veuriation  of  as  well  as  its 
magnitude  very  well. 

A data  run,  taken  for  a F value  of  1.2  is  shown  in  Fig.  All.  As  predicted 
in  Pig.  A5  the  variation  of  Ij.  with  zone  position  is  less  pronounced  for  the 
more  weakly  focussed  case.  Also  shown  are  the  steady  state  vEdue  of 
predicted  by  the  expression  for  which  also  indicates  the  smaller  degree 
of  temporal  variation  expected  at  lower  F values. 

The  resiilts  shown  in  Figs.  AlO  and  All  are  seen  to  be  restricted  to 
values  of  above  about  0.3  with  a corresponding  restriction  on  the  range 
of  the  relevant  thermal  blooming  parameters  studied.  To  expand  this  range, 
measurements  were  taken  with  the  apparatus  kinematics  essentially  fixed  and 
with  the  value  of  P varied  between  5.6  and  l6  W and  at  between  0.21  and 
0.1t5,  €Uid  the  data  from  these  runs  is  shown  as  oscillograms  in  Fig.  A12  and 
in  the  steady  state  parameter  plot  in  Figs.  A13. 

Inspection  of  the  data  shown  in  Fig.  A12  reveals  that  an  increase  in 
the  product  a and  P from  .011  to  0.72  W/cm  does  not  affect  the  ability  of 
the  medium  to  respond  to  the  presence  of  the  laser  beam  and  the  kinematics 
of  the  encounter.  The  relatively  smell  temporal  variation  of  the  peak 
intensity  as  observed  in  Fig.  A12a  is  due  to  the  fact  that  the  df  >arture 
of  Ir  from  unity  is  relatively  small  under  these  conditions.  The  data  in 
Fig.  A12  is  reduced  and  plotted  in  Fig.  A13  as  a function  of  Nl*  this 
plot  only  Ij.  values  occuring  after  the  initial  transient  ere  plotted.  This 
data  is  compared  in  each  figure  with  the  steadj’  state  data  of  Fig.  Al.  In 
each  comparison  the  moving  zone  data  is  seen  to  correlate  with  the  previously 
existing  steady  state  data  within  the  experimented  reproducability  exhibited 
by  the  latter.  As  a result  the  conclusion  is  drawn  that,  for  the  eucovinter 
kinematics  considered  in  the  measurement,  1^  , once  having  reached  a vedue 
corresponding  to  the  steady  state  value  Ij.5s  associated  with  the  thermal 
blooming  parameter  I^q)>  will  cheuage  in  accordance  with  the  sensitiv- 

ity of  Iras  to  the  time  varying  value  of  Nl.  It  is  important  to  emphasize 
that  Ij.  is  determined  by  Nl  only  after  a steady  state  condition  is  reached. 

For  example,  the  data  in  Fig.  A12  correspond  to  a situation  in  which  a 
stagnation  zone  passes  into  the  range  from  behind  the  target.  This  occurs 
at  roughly  0.1  sec  in  the  oscillograms  shown.  Prior  to  that  time,  no 
stagnation  zone  is  present  and  Ij.  is  given  quite  acc\pately  in  terms  of 
Gebhardt  and  Smith's  thermal  blooming  parameter  The  data  shown  in 

Fig.  A13  were  all  taken  at  times  equal  to  or  greater  than  0.25  sec,  leaving 
a time  interval  of  0.15  nec  in  which  Nl  overestimated,  by  as  great  as  an 
order  of  magnitude,  the  severity  of  thermal  blooming.  This  interval  is  seen 


li» 
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from  Fig.  A12  to  be  Independent  of  the  responsivity  of  the  medium,  which 
is  proportional  to  oP,  end  this  indicates  that  the  time  re<iuired  to  set  up 
the  steady  state  condition  is  established  by  the  kinematics  of  the  cell  as 
expressed  by  some  elationship  between  tSz  and  the  ratio  of  the  beam  radius 
ex  some  point  in  the  range  to  w. 

D.  STEADY  STATE  PARAMETERS 

In  the  preceeding  paragraphs  three  distinct  stagnation  zone  parameters , 
designated  Nl'.  Nq  and  N^,  have  been  used  in  analyzing  the  data.  Each  param- 
eter arises  from  a different  approximation  for  the  relationship  at  steady 
state  between  t3z  and  the  ratio  between  the  beam  radius  at  some  point  in  the 
range  to  w.  In  the  case  of  Nt ' the  beam  radius  was  that  associated  with  the 
end  of  the  reuige,  a^-,  (a^  = aTz^)),  where,  in  a fociised  beam,  the  radius  is 
smallest  and  the  contribution  to  thermal  blooming  is  largest.  Although  intu- 
itively Sf  would  thus  seem  xo  be  the  critical  radius  in  determining  the  onset 
of  the  steady  state,  it  wsis  seen  that  becomes  less  effective  in  describing 
steady  state  thermal  blooming  by  a stagnation  zone  as  the  degree  of  focus 
increased.  In  contrast,  is  derived  using  as  the  critical  beam  radius, 
and  Nq  uses  ag , (a^  = (a(z')).  Both  and  Nq  were  seen  tc  exhibit  the  correct 
dependence  on  the  degree  of  focus  with  the  data  exhibiting  a slightly  tighter 

lit  to  an  analytic  function  when  expressed  in  terms  of  N, . 

h 

It  is  interesting  at  this  point  to  try  to  decide  which  of  the  two 
parameters  more  accxirately  describes  the  physics  involved  in  stag- 

nation zone  thermal  blooming.  In  deriving  Nq,  Berger'^'  described  the 
onset  of  siaady  state  situation  with  the  equation 

W(22-2,)tj5«  2(0  U,)d- 0(22))  (A-IV-1) 

where  and  7,2  are  the  boundaries  of  the  ’’thin"  thennal  lens.  Assuming 
a linear  taper  for  the  focused  laser  beam  tho  right  hand  side  of  Eq.  (ATV-l) 
is  equad  to  Uas  and  a pjot  of  the  experimentally  observed  value  of  1^33, 
the  time  elapsed  at  the  onset  of  the  steady  state,  as  a function  of  lies/zfW 
would  be  expected  to  be  a straight  line  whose  slope  is  the  ratio  Zf/(z2-zi). 
Figxure  II4  presents  such  a plot  using  Bergers  data.(^)  Berger  took  his  data 
parame^;lc  in  w at  fixed  z',  and  the  data  points  denoted  with  a cross  represent 
the  run  at  each  z'  for  which  w was  largest.  It  is  seen  that  these  points  fall 
closest  to  the  straight  line  for  which  equals  5.  When  the  slew  rate 

is  large,  transverse  convection  is  most  likely  to  dominate  the  distribution  of 
abaoi'bed  laser  energy  in  the  medium  and  tgg  will  most  likely  reflect  the  rate 
at  idiich  the  slew  generated  transverse  flow  removes  heated  medium  from  the 
beam  path.  At  slower  slew  rates  other  energy  disposal  channels  such  as  bouyant 
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convection  increase  in  in5)ortance,  and  ie  expected  to  fall  below  the  value 
expected  for  the  slew  dominated  case.  As  seen  in  the  figure,  it  would  seem 
that  Berger's  value  of  5 for  Zf/(z2-zi)  is  a reasonable  lower  brand  for  the 
slew  dominated  case.  The  Nq  model  requires  that  the  quantity  tgg  z^w/ag  be 
independent  of  z'  and  have  a value  of  about  20.  However,  when  plotted  against 
z'/z^,  as  in  Fig.  15,  this  quajitity  is  seen  to  show  a systematic  variation, 
increasing  linearly  with  z'/Zj-,  as  ag  decreases.  On  the  other  hand,  when 
plotted  against  z'/z^  as  in  Fig.  l6  for  the  data  of  largest  slew  rate,  the 
quantity  t_„  z,w/a  is  essentially  constant  with  a value  of  about  10  as  an 

5 S 2 O 

upper  bound,  a value  consistent  with  the  data  shown  in  Fig.  17  which  presents 
tsg  vs  U Sq/z^w.  These  results  infer  that  the  vsdue  of  tgg  is  determined  by 
a^  rather  than  by  Sg  and  is  given  by 

4G  (A-IV-2) 


where  the  value  of  f is  roughly  2/5.  In  addition  it  would  seem  that  rather 
than  N-  is  the  relevant  parameter  in  describing  laser  bean  propagation  through 
stagnation  zones. 
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SUMMARY 


In  this  study  measurenents  have  been  made  of  the  reduction  in  peak  on- 
target  intensity  due  to  thermal  blooming  in  stagnation  zones.  Measurements 
were  carried  out  for  both  stationary  and  moving  stagnation  zones.  The  work 
involving  stationary  zones  served  to  expand  the  data  base,  initiated  by 
Berger,  to  Include  a wider  range  of  values  of  F,  at  and  P.  It  was  found  that 
the  stationary  zone  data  was  expressable  in  terms  of  a single  thermal  blooming 
parameter,  N^,  which  is  a slight  modification  on  Berger's  parameter  Nq.  The 
measurements  with  moving  zones  indicated  that  the  measured  intensity  tends  to 
track  the  stationary  zone  value  as  given  in  terms  of  with  the  observed 
ten^joral  dependence  in  the  intensity  being  determined  by  the  temporal  depen- 
dence of  N^. 

One  area  yet  to  be  addressed  in  this  study  of  moving  stagnation  zones 
relates  to  the  question  of  Just  how  fast  are  the  zones  studied  moving  and  how 
fast  are  the  zones  typically  encountered  in  the  naval  scenario  expected  to 
move.  The  quantity  relevant  to  this  question  is  obtained  from  considering  the 
position  of  the  stagnation  zone  within  tbe  range  as  a function  of  time  which 
for  the  laboratory  apparatvis  is  given  by  Eq,  (AIl-l).  Routine  calculations 
reveal  that  dy  /dt  is  largest  for  time  zero  and  when  the  zone  is  Just  enter- 
ing the  range. 


^yp  \ . 

'moximum  "^o 


<AV-1) 


For  the  naval  scenario  considered  in  Ref.  1,  this  value  of  yQ(max)  can  be 
calculated  to  be  O.lt  sec'^.  For  the  data  in  Fig.  A12  this  qiuuntity  is  nine 
aec~l.  Indicating  that  the  zones  studied  are  moving  quite  rapidly  as  viewed 
in  the  context  of  the  naval  scenario.  Thus,  the  conclusions  reached  concerning 
the  stagnation  zones  studied  in  these  experiisents  are  applicable  to  those 
encountered  in  the  naval  scenario. 

In  final  suamazy, conbining  the  results  of  the  present  measurements  with 
previous  experimental’  and  analytical  results  leads  to  the  following  conclu- 
sions: For  a stationary  stagnation  zone,  the  intensity  experiences  an  initial 
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transient,  after  the  laser  is  txirn'^d  on,  before  reaching  a steady  state  value. 
The  tine  elapsed  before  the  latter  is  reached  is  given  approximately  by  10  a^/ 
where  aQ  is  the  beam  radius  at  the  l-aser,  Zf  is  the  laser- to- target  dis- 
tance, and  w is  the  laser  beam  slew  rate.  The  steady  state  intensity,  normal- 
ised to  its  value  in  the  absence  of  thermal  blooming  effects,  is  given  by 

1,  « (l  3.0N|.>f  0.37  N*)  (A-V-II) 

where  the  stagnation  zone  thermal  blooming  parameter  is  given  by 


where  the  symbols  used  stand  for  the  following  physical  quantities: 

= medium's  refractive  index  temperature  coefficient 
o = medium's  absorption  coefficient 

P 5 laser  beam  power 

x'  = distance  of  stagnation  zone  from  the  laser 

a(z)  s laser  beam  radius  « a^  - (a^  - ap)  z/tf 

a^  s bean  radius  at  target  . 

p s medium's  mass  density 

Cp  = medium's  specific  heat 

a^j/aj.  = degree  of  focus. 

For  a moving  stagnation  zone,  the  time  dependence  of  x’  and  therefore  of  is 
determined  by  the  background  wind  speed  and  w,  and  it  was  found  that  the 
above  expression  for  was  obeyed  for  both  stationary  stagnation  xones  and 
zones  moving  at  least  to  ten  times  more  rapidly  than  those  expected  in  a 
typical  naval  scenario. 


(l»2‘/2y)0o 

irpCpuioU'l^O, 


(A-V-III) 
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APPENDIX 

TIME  DEPENDENT  MODEL  FOR  STAGNATION  ZONE  THERMAL  BLOOMING 


A.  INTRODUCTION 

In  section  A-III  and  Ref.  5,  models  were  described  which  account  for 
thermal  blooming  in  the  presence  of  a stagnation  zone  only  after  the  cessation 
of  the  transient  intensity  decay  leading  to  a steady  state.  The  problem  of 
thermal  blooming  in  the  presence  of  a moving  stagnation  zone  is  inherently 
time  dependent,  and  the  development  of  a model  accounting  for  this  time  depen- 
dence is  desirable. 

B.  TIME  DEPENDENT  MODEL 

In  principle  this  need  is  filled  by  a solution  to  Eq.  (A-lII-U)  which  was 
solved  for  the  stationary  zone  by  making  certain  assumptions  about  the  thermal 
lens  thickness  6z.  For  a moving  zone  such  a treatment  is  complicated  by  the 
fact  that  the  boundaries  defining  6z  are  moving  in  a complicated  manner,  and 
although  this  motion  can  be  accounted  for  in  a numerical  solution  to  Eq, 
iA-III-M,  it  was  desired  to  reduce  the  problem  if  possible  to  a simple  one 
parameter  expression  analogous  to  Eq.  (A-III-8)  with  the  elapse  of  time 
included.  Ibis  goal  was  achieved  by  a synthesis  of  a solution  of  the 
thermal  blooming  equation, 

^•(I//A(/7^ds+^l))=-al.  (A-A-i) 

and  some  results  of  the  thin  thennal  lens  analysis  of  thermal  b].ooming. 

Briefly,  using  certain  assumptions  concerning  the  variation  of  I along  the 
propagation  path,  Eq.  (A-A-l)  is  80lved(3)  for  situations  in  which  the  absorbed 
power  density  al  is  balanced  by  a single  energy  loss  channel  for  the  following 
loss  channels:  1)  medium  heating,  2)  thermal  conductivity,  3)  trauisverse 
convection.  For  each  channel  operating  alone,  the  solution  of  Eq.  (A-A-l)  can 
be  expressed  in  the  form 

(A-A-2) 

p - p 

T.<faere  T is  the  product  of  i /a^  , the  fractional  change,  relative  to  ambient 
u^T/u  of  the  index  of  refraction  on  the  laser  beam  axis,  a factor,  f(a£), 
a factor  g^(F)  and,  for  the  case  of  convection,  a factor  ^ which  is  a function 
of  4>  the  ratio  of  the  wind  velocity  at  l to  that  at  the  laser,  'Tq.  These 
quantities  are  listed  in  the  table  below. 
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Loss  Channel  ■ 

AT 

Ua£) 

♦ (<#>) 

Medium  Heating 

aPt  • 

F(F+2) 

irpCpOQ^ 

^2(e-Q^-i  + a/) 

3 

Thermal  Conduction 

OF 

4-rrk  ) 

> a 2/2 

2F(F/nF-(F-|)  ) 
(F-|)2 

Transverse  Convection 

QP 

) 

2FV'i)-/nF) 

2{,pjn<p-{<p-\)) 

Oq^o  I 

1 

(F-I)2 

(^-1)2 

When  all  three  loss  channels  are  included  in  the  energy  balance  equation. 


^ ^ jOCpVAT  . gp  ^ 

a2  t a IT  Q 2 


th-'  value  of  hT  can  be  expressed  as  follows; 


AT* 


aP 


ir^CpQ' 


\t  a2/>cp  o) 


(A-A.3) 


(A-A-1*) 


and  it  is  seen  that  AT  has  a form  similar  to  that  for  the  medium  heating  case 
with  t replaced  by  t*,  a "reduced  time"  given  by  the  quantity  in  brackets  in 
Eq.  (A-A-1*). 


As  a parallel  to  Eq.  (A-A-i*),  it  is  postulated  that  the  general  expression 
for  V in  th'  ‘■agnation  zone  prob?.em  is 


i 


f 

I 


ir^/>Cp  Qq* 


0 mf>CpgcondtP) 


VCmh^F) 

OmgconvCn 


(A-A-5) 


where  tne  factor  (-y^api^  f (o£)e ,g(F)/(iniP^ap^) )"^  is  the  transient  thermal 
blooming  time  as  referred  tc  in  Section  A- II.  It  is  seen  that  ♦(♦)  is  set 
equal  to  2,  its  value  for  a stagr.ation  rone  at  the  ■target  in  Eq.  (A-A-5)*  This 
is  done  because  in  the  presence  of  a stagnation  zone  at  any  other  point  ■within 
the  range  the  transverse  medivm  velocity  ■vanishes  at  the  tone  center  and  the 
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integral  defining  the  value  of  has  no  meaning.  It  is  also  seen  that 
within  the  parentheses  the  beam  radius  is  denoted  by  an',  it  is  concerning  this 
variable  that  the  second  postulate  of  the  model  is  made.  Prom  Eq,  (A-IIT~5) 
it  is  possible  to  determine  the  value  of  z’  at  which  the  placement  of  a 
negative  lens  of  arbitrary  focal  length,  f,  has  the  greatest  perturbative 
effect  on  Ij.,  This  is  conveniently  done  by  determining  the  value  of  z' 
corresponding  to  the  maximum  in  the  second  term  within  the  brackets  in  Eq. 
(A-IlI-5)  and  it  is  found  that 


' 2F-3 

2fn  * • (A-A-6) 

2F-2 

It  is  now  posttalated  that  for  F greater  than  1.5,  set  eq’oal  to  a{zjjj), 

for  F smaller  than  1,5  Sjji  "be  set  equal  to  a^;  in  addition,  the  value  of  v is 
set  equal  to  its  value  that  at  Specifically,  v is  the  product  of  w and 

the  distance  between  z'  and  Zj^.  In  terms  of  the  appeiratus  dimensions  and 
kinematics. 


V 


(A-A-7) 


The  dependence  of  gjQjj(F)t*  on  x and  F is  shown  in  Fig.  AlA.  In  this 
figure  the  value  of  t is  taken  to  be^infinite  so  that  t*  is  composed  only  of 
the  convective  and  conductive  terms.  As  was  the  case  for  the  thermal  lens 
model,  the  effect  of  increasing  F from  1.2  to  5.0  is  to  cause  an  increased 
icj  dependence  of  Y*.  Unlike  the  results  shown  in  Fig.  A-5  for  x(xi»-^»F)t 
further  increase  in  F results  in  a flatter  Xj^ profile.  The  cause  of  this  is 
the  assumption  of  constant  a^  in  the  curves  shown  which  results  in  smaller 
Sjjj  with  increasing  F and  thus  an  increasing  role  of  conduction  rather  than 
convection  in  determining  t*  at  the  larger  F values.  Conduction  was  not 
included  in  the  thermal  lens  model,  and  its  Inclusion  in  the  present  analysis 
points  out  its  applicability  to  situations  other  than  stagnation  zone  thermal 
bloomi.ng  as  shown  by  the  extension  of  the  curves  to  values  of  Xj^  smaller  then 
-i  •‘m/sec  for  which  the  value  of  or  Mq  is  undefined.  Indeed  in  the  absence 
of  & stagnation  zone  and  conduction  and  transient  effects,  the  expression  for 
is  seen  to  be  quite  similar  to  the  parameter  N^S^.used  to  cheuracteri ze 
thermal  blooming  in  the  presence  of  a transverse  wind. 

As  was  the  case  for  N^,,  the  in^ortant  results  of  this  model  is  the 
development  of  a single  parameter  T*  to  characterize  thermal  blooming. 


C.  EXPERIMENTAL  RESULTS 


IRie  sensitivity  of  Ir  to  V*  is  shown  in  Fig.  A2A  for  both  the  present 
experiments,  Berger's  experiments ( 5 ) and  the  analytical  results  of  Ulrich(5)^ 
Wallace(5),  and  Fleck(5»7).  it  is  seen  that  all  the  experimental  results 
correlate  fairly  well  with  y*  and  with  each  other  but  that  the  analytical 
results  for  leurger  P values  fall  away  from  the  lower  F results , Indicating 
that  y*  does  not  conttUn  the  correct  dependence  on  F.  The  moving  zone  results, 
shown  in  Fig.  A-3-A,  indicate  that  the  time  dependence  of  ly  is  determined  by 
that  of  y*  as  was  the  case  for  the  parameters  Nq  and  Nj,. 

The  fact  that  y*  exhibits  incorrect  sensitivity  to  F can  be  altered  by 
simply  replanning  with  b.q  in  Eq.  (A-A-5).  However  consistency  requires 
that  doing  so  forces  Zj^  in  Eq.  (A-A-7)  to  be  replaced  by  zero,  thereby  render- 
ing V and  therefore  y*  time  independent  in  the  limit  of  large  values  of  t, 
also  losing  for  y*  its  sensitivity  to  X|^  and  therefore  stagnation  zone  position. 

In  siaamary,  the  model  described  produces  a thermal  blooming  parameter  y* 
which  correctly  traces  the  temporal  evolution  of  ly  in  the  presence  of  moving 
or  stationary  stagnation  zones  except  for  large  values  of  F.  To  give  the 
correct  dependence  on  F it  would  seem  that  the  third  term  in  psu:enthesi8  in 
Eq.  (A-A-5)  must  be  divided  by  F, 
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APPARATUS  FOR  MOVING  STAGNATION  ZONE  EXPERIMENT 
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FIG.  A~6 
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FIG.  A-3A 
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AEROSOL  EFFECTS  ON  COg  LASER  PROPAGATION 

SECTION  BI 
INTRODUCTION 


Understanding  the  propagation  characteristics  of  high  power  CW  laser 
‘beams  throvigh  an  aerosol  ladened  medlvmt  is  important  for  Na'vy  applications  of 
lasers  for  two  basic  reasons.  First » if  the  laser  is  operating  at  a middle 
infrared  wavelength  for  which  there  is  an  atmospheric  window,  the  effective 
absorption  coefficient  of  suspended  aerosol  can  dominate  over  the  ambient 
molecular  absorption  coefficient  and  second  is  the  problem  of  propagation 
through  fog.  When  these  situations  arise,  thermal  distortion  of  the  high 
power  beam  must  be  calculated  from  a theory  based  on  the  absorption  coeffi- 
cient of  the  aerosols.  It  has  been  demonstrated  experimentally^*^  and 
theoretically^  for  the  steady  state,  CW  case  that  the  propagation  models  of 
molec\ilar  Induced  thermal  blooming  were  applicable  to  the  aerosol  case  by 
merely  replacing  the  gas  absorption  coefficient  with  that  indicative  of 
the  suspended  aerosols.  This  modeling  is  adequate  for  particles  that  cannot 
be  vaporized  by  the  high  power  beam.  The  other  aerosol  propagation  problem,  i.e., 
when  vaporization  occurs,  has  not  been  understood  as  well.  This  second 
important  aerosol  propagation  condition  exists  when  the  target  is  shrouded 
by  fog.  The  detrimental  effects  of  fog  sure  of  concern  whether  the  laser 
device  is  operating  in  the  middle  infrared  region  or  at  longer  wavelengths 
such  as  10.6  microns.  The  thermodynamics  of  the  fog  droplets  have  been 
studied  ^»^and  some  of  the  propagation  effects  have  been  investigated,®*^ 
but  this  work  typically  neglects  the  problem  of  thermal  blooming.  Further- 
more, early  attempts  to  experimentally  measure  fog  induced  thermal  bloomlng^^ 
failed  to  provide  measuiable  distortion,  presumably  due  to  high  transverse 
velocities  that  could  not  be  alleviated.  Therefore,  the  goal  of  this 
research  program  was  envisioned  as  filling  this  information  gap  by  operating 
a high  power  laser  in  a scaled  regime  where  fog  induced  thermal  blooming  could 
be  measured  and  analyzed.  We  have  accomplished  this  goal. 

Another  isgjortant  goal  of  this  research  program  was  to  operate  the  fog 
^.nduced  thermal  blooming  experiments  in  a parametric  arrangement  that  was  a 
scaled  subset  of  conditions  that  might  be  encountered  during  a full-scale 
field  situation.  To  achieve  this  goal  we  have  utilized  a multi-kilowatt, 

CW  CO2  laser  to  provide  intensity  levels  similar  to  those  anticipated  for 
field  conditions.  In  addition,  the  beam  'vras  translated  through  the  fog  cell 
to  simulate  a uniform  wind  transverse  to  the  trajectory.  In  the  propagation 
cell,  a very  concentrated,  artificially-generated  fog  was  used  to  model 
natural  fog  concitions  over  long  atmospheric  paths,  irtiile  maintaining  a droplet 
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size  siailar  to  natural  conditions.  The  parameters  of  the  scaled  experiment 
are  compared  with  approximate  field  conditions  in  Table  B-I.  The  entire  range 
of  parameters  used  during  the  experiments  is  listed  in  Table  B-II. 

This  report  summarizes  our  efforts  to  measure  fog  induced  thermal  distor- 
tion during  the  conditions  discussed  above  (Section  Bill)  and  to  model  the 
physics  of  the  interaction  between  high  power  beam  and  a medium  shrouded 
in  fog  (Section  BII).  Computer  codes  for  thermal  distortion  in  a saturable 
absorber  and  for  Hie  scattering  functions  were  developed  under '^jnited 
Technologies  Corporate  IR&D  fvinds  in  support  this  resea'^ch. 
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SECTION  BII 
THEORETICAL  MODELING 


A.  INTRODUCTION 

Hiis  section  reviews  the  basic  effects  of  a water  droplet  on  propagation 
at  10.6  microns.  The  effects  are  diacussed  in  terms  of  the  extinction  coef- 
ficient, the  absorption  coefficient  and  the  scattering  coefficient.  Pertinent 
examples  of  Mie  calculations  arc  included.  Furthermore,  we  introduce  the 
decay  of  the  water  droplet  as  a consequence  of  high  power  laser  irradiation. 
For  conditions  when  droplet  vaporization  is  significant,  ve  model  the  fog  as 
a saturable  absorber.  The  saturable,  fog  absorber  model  has  subsequently 
been  included  in  our  CW  thermal  blooming  conputer  code.  The  ensuing  results 
of  reduced  thermal  blooming  and  enhanced  power  transmission  are  examined  in 
a limited  sensitivity  study  of  the  code  predictions.  Finally,  with  the  goal 
of  an  eD.pirical  model  in  mind,  we  indicate  an  approach  to  simplified  scaling. 

B,  PROPAGATION  THROUGH  NONVOLATIU:  AEROSOLS 

The  irradiance  loss  on  target  for  linear  propagation  is  given  by 


I 

JS* 


a,^,U)dz 


(B-1) 


where  L is  the  pathlength  through  a fog  of  extinction  coefficient  It  is 

assumed  that  the  beam  exerts  no  modifying  influence  on  the  fog  droplets.  If 
the  fog  is  uniform  and  homogeneous  then  one  may  write  Eq,  (B-l)  sioply  as 


(B-2) 


The  extinction  coefficient  is  given  by 


f<0 

»r*0„^(f)N(r)dr 

0 


(B-3) 
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in  which  is  the  extinction  efficiency  for  droplets  of  radius  r with  a 

sise  distribution  N(r),  For  a aonodlsperse  droplet  population  Eq,  (B-3)  is 


a,*t«7rr*0„t(nc 


(B-I4) 


where  C is  the  concentration  of  fog  droplets  per  unit  volume. 

In  this  development  it  is  important  to  identify  two  mechanisms  responsible 
for  intensity  loss  during  the  aerosol  propagation  problem.  These  are 
absorption  and  scattering,  such  that 


Qwt(»’)«0ob»(0-^0,co  (r)  (B.5) 

The  scattering  coefficient  has  no  bearing  on  thermal  distortion  considerations 
except  with  respect  to  a general  intensity  loss  as  the  beam  propagates  to  the 
target.  However,  the  absorption  coefficient  is  directly  responsible  for 
transferring  energy  from  the  laser  beam  to  heating  of  the  ambient  gas.  These 
induced  temperature  chauiges  in  the  gas  can  cause  severe  beam  distortion  and 
are  the  result  of  nonlinear  propagation  effects. 

Now,  if  one  is  concernel  with  nonlinear  propagation  effects  and  assumes 
that  the  droplets  remain  unmodified  by  the  propagating  beam,  the  amount  of 
thermal  blooming  is  conveniently  modeled  by  the  peak  target  irradiance  rela- 
tive to  the  undistorted  value.  The  empirical  expression^  for  the  reduced 
target  irradiamce  (considering  only  thermal  distortion)  is 


where 


(B-6) 


Nf« 


-'nr  aqbi  PL*  , 


(B-7) 
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The  basic  parameters  are  defined  as 


tjt  “ index  change  per  unit  terperature 

p , i 

P = laser  power  v 

n = unperturbed  index  of  refraction 
p » aablent  gas  density 
Cp  = ambient  gas  heat  capacity 
V = transverse  wind  velocity 
a * 1/e  intensity  radius 


In  Eq.  (B-7)  ^(®^ab  correction  factors  for  large  absorption 

coefficients  (and/or  long  path  lengths)  end  strong  focussing: 


(B-8) 


(B-9) 


F is  the  Fresnel  nianber.  Note  in  this  development  that  the  absorption  coef- 
ficient (assuming  negligible  background  molecular  absorption)  must  be  repre- 
sentative of  the  fog  droplets. 


aob.*’^«’*Oobs('r)C 


(B-IO) 


Therefore,  a fundenental  part  of  the  linear/nonlinear  fog  propagation  problem 
consists  of  identiiying  the  absorption  and  scattering  efficiencies  of  t^ie 
droplet  'Constituents  at  the  wavelength  of  interest. 
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A Mie  scattering  computer  program  has  been  developed  at  UTRC  and  was 
utilized  to  provide  the  necessary  efficiency  factors.  This  program  calculates 
the  angular  scattering  function  for  a droplet  of  given  size.  Two  examples  are 
given  in  Fig.  B-1,  These  curves  are  for  polarization  perpendicular  to  the 
plane  of  incidence.  In  the  case  of  the  0,5  wm  radius  droplet  we  observe 
nearly  perfect  Rayleigh  scattering^^  while  in  the  5 pm  radius  case  the 
angular  function  has  developed  more  structure.  This  angular  scattering  func- 
tion is  integrated  over  all  angles  to  derive  the  total  scattering  efficiency' 
for  a droplet  radius.  The  absorption  efficiency  is  also  calculated.  Figure  B-2 
gives  the  ratio  of  the  absorption  and  scattering  efficiencies  for  droplet  sizes 
in  the  range  of  interest  for  this  project.  Furthermore,  the  absolute 
absorption  efficiency  is  illustrated  in  Fig.  B-3,  The  calculations  illus- 
trated in  Figs.  B-2  and  B-3  can  be  used  to  estimate  the  severity  of  the  linear 
and  nonlinear  propagation  problems  through  a nonvolatile  fog. 

The  parameters  for  a full-scale  propagation  example  in  light  fog^^  are 
listed  in  Table  B-III.  The  target  irradiance  losses  are  substantial  with  the 
effect  of  fog  induced  thermal  blooming  being  most  severe.  In  contrast,  only 
nominal  losses  are  suffered  for  the  clear  atmosphere.  In  summaiy,  this 
example  alone  serves  as  suffici,  .t  Justification  for  concern  over  the  degrading 
effects  that  occur  during  foggy  conditions. 

Note  that  this  discussion  was  based  on  the  assumption  that  the  fog  was  not 
modified  by  the  laser  beam.  In  reality  a fog  droplet  is  very  volatile  and 
under  certain  conditions  a sufficiently  powerful  laser  may  destroy  a large 
portion  of  the  droplet  pop’ilation.  In  this  situation  the  outlook  for  high 
power  propagation  is  not  ntiariy  as  bleak.  The  theoretical  results  for  prop- 
agation in  actively  modified  fog  are  discussed  in  the  next  section. 

C.  PROPAGATION  THROUGH  ACTIVELY  MODIFIED  FOG 

Figure  b-1*  presents  the  calcxxlated  decay  of  a 10  pm  fog  droplet  irradi- 
ated by  a high  power  beam.  In  this  illustration  the  droplet  enters  a uniform 
5 kW  beam  with  a 1.5  cm  radius.  A scaled  wind  of  15  cm/s  is  blowing  trans- 
verse to  the  trajectory.  By  the  time  the  droplet  has  reached  the  beam  axis 
it  has  disappeared.  The  Jfl5)crtance  of  the  droplet  lifetime  problem  is  com- 
pounded when  Eqs.  (B-l*)  and  (B-IO)  are  recalled  In  which  the  extinction  and 
absorption  coefficients  are  a strong  f\in.'!tion  of  the  fog  droplet,  size.  As 
a consequence,  ^xjdels  for  power  transmission  and  thermal  distortion  must  be 
appropriately  modi  fled  to  reflect  the  influence  of  the  decaying  droplets.  In 
this  section  we  address  such  a roodiflcatioin  with  an  initial  discussion  of  the 
droplet  thermodynajnics. 
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Many  references 


6,9,15 


give  the  rate  equation  for  a vaporizing  droplet 


d(4/3’rr*^Hv) 

dt 


-Airv 


lOobt 


(B.-il) 


In  this  development  it  is  assumed  that  the  droplets  are  not  vaporized  fast 
enoxigh  to  force  the  velocity  of  the  escaping  vapor  to  exceed  the  acoustic 
velocity®.  In  addition,  Mullaney,  et  al9  show  that  the  solution  to  this 
equation  may  be  written  as 


rCt.Msfj# 


At/T® 


(B-12) 


where 


, , KqW^tV) 


(B-13) 


for 


o ■ bulk  absorption  eoeffl  :lent  of  H-0 
■ thermal  conductivity  of  air 
R * ijas  constant 

T(?)  r temperature  et  drop  surface 
* HgO  heat  of  vaporization 
D ^ diffusivity  of  HgO  vapor  in  air 
w « molecular  wei^t  of  H^O 
P^fr)  ■ partial  pressure  of  water  vapor  at  !r(r) 
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From  Rer.  6 we  find  that  for  conditions  of  interest  T(r)  s ambient,  an 
approximation  good  to  abo\it  10  percent.  Hence,  P^(r)  is  the  vapor  pressure 
of  the  saturated  solution  at  the  ambient  temperature ..  For  the  following  con- 
stant values 


6 


X 10^  cm  ^ 


0*1  g/cm^ 

H “ 2U48  J/gm 

V 

-k 

K » 2,55  X 10  J/csnSecoC 

a 

R * 8.318  J/°K  mole 
.0 


T = 293  K 
C 


-1  2 

2.k  X 10  cm  /sec 


rf  = 18  gm/mole 


wv 


» 23.1*  Bb 


tI  3 11  J/cm  with  I in  w/cm  , so  we  rewrite  £q.  (B-12)  as 


rd.C'^rjff 


rAti/  II 


(B-lU) 


In  general,  we  e:q)ect  Eq.  (B-Jl*)  to  be  a reasonable  solution  to  Eq.  (B-ll)  for 
droplets  with  radii  less  than  about  10  pn. 

Furthermore,  ve  have  performed  a linear  refression  analysis  of  os  a 
function  r(cm)  and  find 


®obs  t*")  • *0“^ 


(E-15) 
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Hie  correlation  coefficient  for  this  fit  to  xhe  Mie  calculated  values  irs 
better  than  99  percent  for  r ■ 10“^  to  1G"3  em,  we  may  further  simplify  this 
expression  (error  less  than  10  percent  in  range  cf  inxerest)  to 


Qobt(f^)=  (5.18  X 10*)  r 


(B-16) 


C<»ibining  Eqs,  {B<-16),  (B-iU)  and  (B-lOi 


3 - »An/  u 


(B-ir) 


It  is  then  obvious  that  a will  depend  on  a droplets  initial  size  and  the  total 
time  it  has  been  iUuinlnated  by  the  high  power  beaa:,  Since  I uiay  be  a func- 
tion of  position,  and  the  meditun  is  moving  with  velocity  v,  ve  may  rewrite 
Eq.  (B-17)  as 


(K)  * (LW X 10*)  irCr.*  e 


(I)-16) 


upon  considering  on^  the  plane  of  the  beam  containizig  the  wind  vector.  Now, 
Gebhardt  end  Smlth^®  show  the  perturbation  in  the  ambient  gas  temperature  is 


ATW* 


I 

pCpy 


(B-l?) 


At  this  point  ve  must  acknowledge  Uiat  a signifir^ant  portion  of  tJie  eiiei’tfir  is 
lost  to  the  vaporization  mochanisu  end  cannot  be  co.;auct«v^  i'c  th*  gas.  Hist 

is 
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[Tctr 


(B-20) 


where  c is  the  fraction  of  absorbed  heat  that  is  conducted  tc  ttiC  gas  from  the 
droplet,  e is  Just^ 


€ • 


KoR*T»(r) 

OH^w^P^yCr) 


(B-21) 


which  has  a value  between  0.3  and  O.U.  We  will  use  0.33.  Finally 


^Cpv 


dx' 


(B-22) 


where  T{s)  is  the  one~diinensional  perturbation  tei!?)erature  of  the  aobient  gas 
when  the  fog  is  subject  to  vaporization.  If  we  nomallze  tliis  teniperaturc 
change  by  that  experienced  for  no  induced  droplet  vaporization  Eq.  (B-22) 
becoiaes 


AT  /-<p— "» 

^ /“  lUldV 

y-OD 


(B-23) 
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■nie  importance  of  Eqs.  (B-22)  and  (B-23)  can  be  illustrated  by  considering 
a uniforjn  transmitter  illumination  function.  Equation  B-22  becomes 


. . (9.05X10®)  wCr?  f 

AT(x)= = J 


{B-2M 


where  I is  the  laser  intensity  and  x is  the  position  in  the  beam  cross  section 
along  t§e  direction  of  the  wind.  Purthe^Tnore , Eq.  (B-23)  becomes 


AX 


(B-25) 


The  above  equations  indicate  that  the  maximum  ten^erature  perturbations  may  be 
severely  limited  due  to  the  saturable  absorber.  The  significance  of  this 
effect  is  displayed  in  Fig.  5 where  we  have  graphed  the  relative  ambient 
temperature  rise  for  a nonsaturable  and  a saturable  absorber.  This  exan5)le 
shows  that  the  on-axis  temperatiire  perturbation  for  the  saturable  fog  absorber 
is  only  about  lU  percent  of  that  for  a nonvolatile  absorber  (or  molecular 
absorber).  Even  more  Inqportant  is  the  fact  that  the  perturbed  temperature 
gradients  are  also  drastically  limited.  This  condition  implies  reduced 
thermal  distortion  of  the  laser  beams  spatial  intensity  distribution. 

The  expressions  corresponding  to  Eqs.  (B-2U)  and  (B-2^)  for  a Gaussian 
illumination  function  are 


AT(x)  * 


(0.33)(l.i8xtO®)Trrfio  r * — " (ift’)  [' 

® 


pCpv 


■dx' 


(B-26) 
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and 


[u  art  (Vo)] 

dx 

2 1 

1 -f  erf(x^)] 

(B-27) 


where  I is  the  peak  intensity  for  a beam  with  e radius  of  a. 
o 

In  conclusion,  these  modified  ten^perature  perturbations  can  be  used  in  a 
thermal  blooming  code  to  calculate  the  resulting  distortion  of  a high  power 
beam  and  to  study  the  importance  of  the  saturable  absorber. 

D.  PROPAGATION  CODE  WITH  SATURABLE  FOG  ABSORBER 

The  saturable  fog  absorber  model  has  been  incorporated  into  a CW  thermal 
blooming  code^*^.  Tlie  modified  code  takes  the  initial  fog  droplet  size  and 
concentration  and  uses  Mie  calculations  tc  develop  the  appropriate  extinction 
and  absorption  coefficients.  The  beam  then  propagates  through  one  incremental 
step  toward  the  target,  forcing  droplet  evaporation.  The  phase  changes  due  to 
index  perturbations  from  the  fog  induced  temperature  perturbations  are  iden- 
tified and  utilized  as  initial  infomation  for  the  next  Z step.  The  back- 
ground molecular  absorption  coefficient  is  eJ-SO  accounted  for.  In  addition, 
the  extinction  from  the  modified  fog  is  calculated  such  that  the  initial 
Intensity  for  the  subsequent  Z step  is  appropriately  modified.  Calculations 
over  the  beam  cross-section  are  performed  on  a rectangular  mesh  that  is 
adapted  to  the  convergence  of  the  focused  beam,  and  the  Incremental  Z stepping 
toward  the  target  is  activcl;/'  adjusted  to  keep  the  accumulated  distortion  per 
step  very  small.  The  general  calcul ational  procedures  are  similar  to  other 
nonlineeur  propagation  codes^^^. 

An  exas^le  of  the  output  graphics  from  the  code  is  given  in  Fig.  B-6. 
Shown  for  typical  experimental  conditions  are  the  phase  distribution,  inten- 
sity distribution,  transverse  intensity  profile  and  "ihe  integrated  intensity 
in  the  target  plane.  This  particular  example  is  for  0.^  kV,  and  thermal 
distortion  of  the  intensity  distribution  is  evidenced  as  is  a large  total 
power  loss  via  extinction.  When  the  power  is  Increased  to  U.O  kW,  Fig.  B-7 
shows  that  thermal  blooming  decreases  and  the  overall  relative  transmission 
to  the  target  plane  has  increased  aubstantially.  In  an  attempt  to  under- 
stand the  scaling  of  these  effects,  we  have  looked  at  some  code  runs  in  a 
more  comprehensive  fashion. 
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To  date  the  new  propagation  code  heis  ’been  exercised  over  a limited 
range  of  parameters.  For  instance,  we  have  performed  a sensitivity  study 
study  of  parametric  variations  on  the  nominal  thermal  blooming  conditions  of 
one  case  from  the  experimental  data  runs  to  be  discussed  in  Section  Bill.  In 
this  study  we  have  investigated  I^el  Eq.  (B-6))  as  a function  of  the 

distortion  parameter  Nf,  when  Hf  is  determined  by  the  absorption  coefficient 
of  the  unmodified  fog.  The  result  is  presented  in  Fig.  B-8.  Table  B-IV 
reviews  the  range  of  parameter  variations.  As  P and  v are  adjusted  providing 
different  "burn-through"  capabilities,  the  breakdown  of  the  classical  thermal 
blooming  model  Eq.  (B-6)  becomes  very  apparent.  Additionally,  the  value  of 
N-  shows  a sensitivity  to  variations  in  the  droplet  radius.  This  is  due  to 
the  strong  dependence  of  the  absorption  coefficient  on  r (Eq.  (B-IO). 

However,  I^gj^  does  not  appear  to  be  very  sensitive  to  the  small  adjustments 
of  r.  Both  and  Nf  show  a much  weaker  dependence  on  the  fog  concentra- 

tion. Ihe  resets  of  a small  variation  in  C reside  between  the  indicated 
points  (I  and  along  the  "P- variable"  line. 

The  enhancement  of  high  power  propagation  through  the  saturable  fog  (and 
breakdown  of  the  classical  model)  can  be  explained  as  occurring  from  two 
effects.  First  of  all,  the  satiurable  absorber  permits  only  limited  tempera- 
ture pertiirbations  to  be  produced  in  the  ambient  gas.  This  results  in  reduced 
spatial  distortion  of  the  bean.  Ibe  second  enhancement  mechanism  is  singily 
that  of  droplet  destruction;  consequently  the  steady  state  effective  absorp- 
tion coefficient  is  significantly  reduced.  As  a result,  overall  transmission 
increeises  in  addition  to  the  inhibition  of  spatial  distortion.  In  the  data, 
these  effects  may  be  separated  by  analyzing  the  target  intensity  from  two 
additional  points  of  view.  We  suggest  that  the  effects  of  reduced  spatial 
distortion  may  be  understood  by  Investigating  the  area  of  the  bloomed  beam 
normalized  by  that  of  the  undistorted  beam  at  the  1/e  intensity  contour. 
Furthermore,  the  mechanism  of  enhanced  transmission  may  be  analyzed  by 
studying  the  total  power  delivered  to  the  target  plane  through  the  1/e  inten- 
sity contoiiT  for  the  distorted  case  compared  to  the  undlstorted  case.  These 
pareuneters  will  be  defined  as  A^^j^  and 

The  Ap^^  results  are  summarized  in  Fig.  B-9.  In  this  exBiig>le  we  have 
maintained  Nf  as  the  abscissa,  but  do  so  recognizing  that  Nf  may  not  be  a 
good  par'ameter  for  the  scaling  of  actively  modified  fog  results.  The  cal- 
culations show  that  as  P increases  or  v decreases  the  fog  is  more  effectively 
destroyed,  and  the  amount  of  spatied  beam  distortion  is  dramatically  reduced 
relative  to  the  classiced  model.  On  the  other  hand,  if  P and  v are  fixed 
and  the  droplet  radius  or  the  droplet  concentration  is  adjusted,  then  the 
amount  of  thermally  induced  distortion  seems  to  have  the  same  slope  or  func- 
tional dependence  on  as  the  gas  model  but  with  a significantly  lower 
magnitude. 


Figure  B-10  illustrates  P^ei  ^ Function  of  the  initially  unmodified 
fog  Large  increases  in  transmission  are  predicted  for  either  increased 

power  or  decreased  wind  velocity.  Also,  Eq,  (B-2)  seems  to  still  he  valid  in 
the  modified  fog  examples  for  fixed  P and  v with  variable  r or  C,  but  a new 
effective  extinction  coefficient  must  be  identified.  In  essence,  this  implies 
that  a 10.6  ym  transmission  through  the  modified  fog  remains  a relatively 
simple  function  of  the  total  remaining  liquid  water  content  in  the  trajectory^^. 
These  code  data  suggest  that  some  siniplified  scaling  may  be  possible  for  the 
nonlinear,  fog  propagation  problem. 

E.  SCALING  AND  A MODIFIED  EMPIRICAL  MODEL 

A reasonable  engineering  approximation  to  the  code  results  illustrated  in 
Fig.  B-10  may  be  developed  by  realizing  that  for  the  droplet  size  rang*  of 
interest  the  extinction  coefficient  is  approxi.metely  equal  to  the  absorption 
coefficient.  With  this  assumption  we  may  utilize  Eq.  (B-17)  to  write 
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(B-28) 


which  becomes 
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For  slB^pliclty,  ve  will  choose  to  calculate  the  scaling  at  beam  center  where 
X • 0 and  include  a fitting  parameter  D such  that 
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(B-30) 


We  select  a value  of  D * 0,U37  to  provide  a good  empirical  fit  between  the 
code  data  and  the  simplified  formulation  of  Eq.  (B-30)  (one  should  note  that 
this  fitting  constant  will  change  depending  on  the  illumination  function  and 
the  Fresnel  number  of  the  beam).  Equation  (B-30)  becomes 
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When  this  modified  extinction  coefficient  is  used  in 


PrtI  ‘ ^ 


(B-32) 


the  sinplifled  model  agrees  to  better  than  10  percent  with  all  but  one  of  our 
code  sensitivity  results.  At  one  low  power  (0.5  KW)  case  where  bum- through 
was  not  veri!^  effective  the  agreement  dropped  to  17  percent. 

Now,  since  s we  may  estimate  a new  effective  absorption 
coefficient  to  be  used  in  calculating  and  of  the  enpirical  model. 

In  other  words 
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-(O.OSS*) 

“thtfmof  * ®ab»  * (B-33) 

blooming  initial 


The  resiilts  of  Fig.  B-3.1  occur  when  Eq.  (B-33)  is  used  to  calculate  a new  Nf 
for  the  conputed  values  of  Fig.  B-9,  TTiis  evidence  iinplies  that  the 

parametric  modeling  by  Eqs.  (B-6)  and  (B-7)  remains  valid  and  that  if  the 
appropriate  modified  absorption  coefficient  is  known  within  the  beafflt  then  the 
thermal  distortion  occurs  as  derived  for  the  molecular  absorption  theories. 

The  one  point  where  our  simplified  scaling  has  not  worked  well  was  for  large 
P/va  (>  100).  In  this  case  the  simplified  model  predicted  a lower  than 

calculated  by  the  computer  code.  Our  work  indicates  that  during  these  con- 
ditions of  large  P/va  the  simple  exponential  law  of  Eq.  (B-28)  must  be  gener- 
alized to  include  an  integration  edong  the  trajectory.  In  such  a case 
scaling  will  become  more  conqpli rated. 

P.  LIMITATIONS  OF  THE  SIMPLIFIED  MODEL 

The  analysis  of  the  computer  calculations  discussed  in  this  section 
represents  only  a limited  sensitivity  study,  and  therefore  does  not  Investigate 
all  limits  of  the  simplified  scaling,  e.g.,  very  large  P/va.  Furthermore,  the 
illumination  function  was  assumed  to  be  Gaussian  with  a fixed  Fresnel  number. 

A change  of  either  of  these  geometric  parameters  will  affect  a change  in  the 
constants  of  the  scaling  equations.  Finally,  the  assumptions  that  X ■ 10.6  ym 
and  r « 2 - 3 um  simplified  the  calculation  of  Q^bs  ^xt* 
deviations  from  these  values,  we  ejqpect  simileur  parametric  effects  to  occ\ir 
but  with  ilfferent  constants  determining  the  amounts  of  burn-through  and 
reduced  tl.ermal  distortion.  In  conclusion,  it  appears  that  the  siinplified 
model  cou.d  be  used  to  predict  high  power  propagation  effects  during  realis- 
tic scenario  conditions  if  X ■ 10.6  ym,  P < 10  and  P/va  ^ 100. 
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SECTION  Bill 

EXPERIMENTAL  MEASUREMENTS 


We  experimentally  investigated  the  possihilltles  of  reduced  thermal 
distortion  and  enhanced  transmission  for  high  power  laser  propagation  through 
fog.  This  portion  of  the  report  summarizes  these  measurements  and  the  equip- 
ment utilized  to  perform  the  study. 

A.  APPARATUS 

1.  System  Overview 

A simplified  schematic  diagram  of  the  expermental  system  is  provided  in 
Fig.  B-12.  The  arrangement  consists  of  six  basic  components:  laser,  uniform 
translator,  cell,  fog  generators,  fog  diagnostics  and  target  plane  diagnostics. 

In  order  to  perform  the  experiment  providing  full-scale  intensity  levels 
it  was  necessary  to  use  a relatively  high  powered  laser.  We  employed  an 
electric  discharge,  transverse  flow  device  developed  and  maintained  at  UTRC 
for  the  U.  S.  Navy. Figvire  B-13  is  a photograph  of  the  laser..  In  the  MOPA 
configuration,  this  device  has  provided  more  than  25  hW;  however,  our  investi- 
gation typically  required  less  than  10  kW  average  output  power.  During  our 
experiments,  the  200-300  watt  oscillator  beam  was  chopped  at  500  Hz  with  a 
50  percent  duty  cycle  to  allow  synchronized  detection  of  fog  propagation 
phenomena.  The  chopping  effect  is  considered  negligible  in  this  experiment 
because  the  fog  is  irradiated  for  sufficient  time  duriu  each  cycle  to  allow 
the  droplet  thermodynamics  to  approach  a steady-state,^  and  the  amount  of 
effective  pulse  overlap  is  large  to  ensure  CW  thermal  distortion  effects^^ 
(pulses  per  flow  time  of  33-113). 

From  the  laser,  the  beam  was  relayed  to  a translating  mirror  via  a tele- 
scope. The  moving  mirror  provided  an  apparent  uniform  wind  vdlth  velocities 
adjustable  between  10  and  50  cm/s.  The  translation  rate  was  calibrated  with 
a digital  system  that  utilized  a stationary  photodiode  viewing  an  LED  through 
a photographic  negative  ettaohed  to  xhe  moving  mirror.  The  negative  had  a 
well  defined  cyclic  bar  pattern  to  generate  regular  signals  in  the  diode  with 
a frequency  dependent  on  the  mirror  translation  rate.  During  calibration 
procedures,  the  velocity  of  the  mirror  was  found  to  be  uniform  to  better  than 
10  percent  over  the  range  given  above.  The  beam  was  then  deflected  into  the 
propagation  cell  by  the  translating  mirror- 
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2.  A«ro8ol  Propagation  Cell 

configuration  of  the  aerosol  propagation  cell  is  shown  in  Fig. 

B-lU.  Ihe  assoahled  chaniber  is  approximately  10-m  long  and  consists  of  five 
equal  length  sections  of  approximately  30>cm-ID  standard  wall  thickness 
aluminum  pipe.  The  recirculation  loop  of  the  test  chamber  consists  of 
five  equal  sections  of  approximately  7.5-cm-ID  aluminum  Schedxd.e  ^*0  pipe. 

All  sections  of  pipe  are  Interconnected  and  sealed  by  means  of  Vitaulic 
couplings.  Each  section  of  the  test  chamber  and  recirculation  loop  is  mounted 
on  its  own  movable  stand. 

Two  fans  for  maintaining  circulation  through  the  cell  are  provided.  A 
7.5-em  dia  529  L/mln  capacity  fan  is  located  in  the  center  section  of  the 
recirculation  loop.  A second  fan,  17.5-cm  dia  and  2175  L/min,  is  positioned 
near  the  cell  input.  These  fans  were  vital  elements  for  stabilizing  the 
relative  humidity  for  fog  absorption  data  runs  and  for  mixing  CO^  with  N-  for 
molecular  absorption  data  runs.  Each  end  of  the  cell  was  sealed  with  a 0.0005 
mil  thick  mylar  window  until  the  appropriate  interior  conditions  were  generated. 
Then  the  windows  were  cut  away  and  a test  conducted.  A pressure  relief  valve 
set  at  approximately  1.007  atm  is  located  on  the  test  chamber  and  prevents 
overpressurization  of  the  mylar  windows  during  introduction  of  the  fog  and 
carrier  gas  into  the  test  chamber. 

Four  sets  of  two  diametrically  opposed  5-cm  dia  viewports  are  provided 
on  the  first,  third,  fourth  and  fifth  test  chamber  sections,  respectively. 

These  ports  enabled  us  to  operate  fog  uniformity  tests  along  the  length  of 
the  chamber • 

Aerosol  generators  are  located  on  all  sections  of  the  propagation 
cell.  Furthermore,  the  second  section  contains  the  scattering  and  extinction 
diagnostics  for  characterization  of  the  fog.  Details  of  the  second  section 
are  provided  schematically  in  Fig.  B-15.  The  third  section  contains  additional 
diagnostics,  e.g.,  spectrophone  head,  relative  humidity  probe  and  temperature 
probe.  These  diagnostics  are  described  in  more  detail  in  following  sections 
of  this  report.  Fig\ure  B-l6  provides  a photograph  of  the  entire  cell,  viewing 
from  the  target  toward  the  input  end. 

3.  Aerosol  Generator 

In  selecting  a droplet  generator  for  this  experiment,  we  kept  four 
goals  In  mind.  First,  it  is  desirable  to  conduct  the  propogation  experiment 
with  a monodlsperse  population  of  droplets.  This  is  important  from  the 
BKSdellng  point  of  view  because  the  propagation  code  as  we  have  developed 
it  to  this  point,  assumes  a monodispersed  population.  Since  Eq.  B-l8  and 
Figs.  B-8  - 10  show  a sensitivity  to  the  aerosol  radius,  we  place  the  impor- 
tance of  a monodi8pe'*se  distribution  as  a primary  goal  in  selecting  a fog 
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generator.  As  a second  goal  we  would  like  to  operate  with  a d-^oplet  radius 
that  closely  resembles  net^-iral  fog  droplecs.  However,  this  goal  must  be 
teu5)ered  with  the  fact  that  oiir  experiments  are  conducted  over  a very 
short  range  compared  to  that,  of  a full-scale  scenario,  and  we  mast  maximize 
the  droplet  absorption  efficiency  in  order  to  transfer  sufficient  energy  to 
the  ambient  gas  to  simulate  realistic  thermal  blooming  conditions.  This  means 
that  the  droplet  size  should  be  the  leurgest  we  cen  generate  compatible  with 
the  natural  scenario  (-^ee  Fig.  T -3).  On  the  other  hand.  Fig.  B-2  indicates 
that  for  large  droplets  a significant  fraction  of  the  intercepted  energy  will 
be  scattered  in  addition  to  being  absorbed.  This  portion  of  the  lost  energy 
is  difficult  to  document  during  the  experiments,  so  as  a practical  matter 
we  would  also  like  to  maintain  as  high  an  eosorption  to  scatterlnp,  ratio 
a:;  poseible.  Droplet  radii  hetveen  2 and  3 microns  represent  a good  compro- 
mise between  these  various  llmivations.  The  third  goal  envisiorAed  as  a 
specification  of  the  aerosol  generator  was  that  it  must  produce  a very  large 
number  of  droplets  per  second.  If  we  require  on  the  order  of  10^  droplets 
per  cm^  and  the  voliome  of  our  chamber  is  7.1  x 10^  cm^,  then  we  must  provide 
about  10^^  droplets  per  data  run.  In  actuality  the  nianber  must  exceed  this 
value  when  taking  into  account  those  droplets  that  will  be  impacted  on  the 
sides  of  the  chamber  during  generation  or  that  will  settle  out^^  Just  prior 
to  a data  inin.  Hence,  to  keep  the  initial  set-up  time  of  the  propagation 
conditions  to  a reasonable  minimum,  the  aerosol  generator  must  produce  about 
lO*^  droplets  per  second.  As  a final  goal  in  our  selection  process  the  fog 
generator  should  operate  at  relatively  low  air  and  water  pressures  so  that 
specialized  pumping  systems  are  not  required. 

Table  V summarizes  the  decision  making  process  in  which  we  were  forced 
to  compromise  the  monodisperse  specification  in  order  to  obtain  a droplet 
generator  satisfactory  for  fulfillment  of  this  investigation's  overall  goal 
of  measuring  and  analyzing  fog  Induced  thermal  distortion.  Of  two  companies 
which  produce  aerosol  generators  which  utilize  the  shearing  forces  of  an 
acoustic  field  to  aid  the  production  of  unifoim  droplets,  we  purchased  genera- 
tors from  the  Sonic  Development  Corporation.  This  coiiq>any  claimed  slightly 
better  adjustability  of  the  droplet  radius  than  the  other.  A photograph  of 
the  generator  is  given  in  Fig.  B-l6. 

After  receipt  of  the  generators,  we  attempted  to  document  the  population 
distribution  emitted  at  the  nozzle  outlet.  Initially,  we  used  a Corporate- 
owned,  pulsed-ruby  holographic  system  to  make  records  of  the  generated  droplets. 
This  technique  was  abandoned  after  a few  etteiqptB  because  the  holographic 
s^'sten  could  only  resolve  droplets  down  to  10  um  in  radius,  leaving  most 
of  the  :populatjon  inadequately  measured.  Then,  we  employed  two  angle 
scattering  in  a second  attempt  to  characterize  the  droplets.  However, 
concentration  fluctuations  in  the  immediate  vicinity  of  the  nozzle  made 
these  muasurements  somewhat  unreliable.  Our  final  effort  to  characterize 
the  fog  produced  by  the  generators  was  conducted  by  filling  a large  chamber 
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(2  IB  long  euQd  0.25  a diaineter)  with  the  fog  in  a aituation  similar  to  that 
of  the  real  experiments.  This  procedure  was  successful,  proved  that  the 
generator  was  adequate  and  is  discxissed  in  the  next  section. 

1*.  Fog  Cheu’acterizatlon 

Since  water  droplets  are  rather  volatile  it  is  desirable  to  have  real 
time  ‘iocumentation,  of  the  mean  aerosol  diameter,  the  fog  droplet  concentra- 
tion, the  extinction  and  absorption  coefficients  at  10.6  um,  and  the  visible 
extinction  coefficient.  As  indicated  by  Ref.  11  these  measurements  may  be 
achieved  with  optical  techniques  that  depend  or.  the  relative  sizes  of  the 
di'oplets  and  the  probe  laser  wavelengths.  We  have  accomplished  this  goal 
of  fog  characterization  by  utilizing  transmission  and  scattering  measurements 
at  10.6  um  and  6328  2. 

Data  from  the  two  wavelength  transmission  diagnostics  provide  most  of  the 
basic  information.  As  an  example,  the  extinction  coefficient  measured  at 
6328  S provides  an  estimate  of  the  visibility  in  the  experimental  fog  con- 
ditions. Reference  23  shows  that  the  visibility  through  a fog  is  simply 
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Using  this  ps^amster  we  may  compare  our  scaled  measurements  with  full  scale 
scenarios  in  tdjich  fog  conditions  are  described  in  terms  of  the  visibility. 
In  addition,  the  transmission  diagnostic  provides  us  with  the  extinction 
coefficient  of  the  fog  at  10.6  wm  which,  as  we  discussed  earlier  in  Section 
BII,  is  a fundamental  parameter  in  our  modeling  of  nonlinear  propagation 
through  fog. 


Furthermore,  Pig.  B-17  shows  that  we  may  measure  the  average  fog  droplet 
radius  simply  by  taking  the  ratio  of  the  extinction  coefficients  at  the  two 
wavelengths  of  our  transmission  diagnostics.  With  this  estimate  of  the 
droplet  radius  we  may  go  bask  to  the  extinction  measurement  at  6328  S and 
calculate  the  droplet  concentration  by  observing 
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Reference  2k  explains  that  this  method  of  measuring  the  aerosol  concentration 
if  reliable  to  about  10  percent.  Finally,  since  we  have  a measure  of  the 
droplet  radius.  Fig.  B-2  can  be  consulted  to  give  an  indication  of  the  absorp- 
tion coefficient  at  10.6  ym  for  use  in  the  theoretical  predictions  from  the 
fog  propagation  models. 
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Wt'  recognize  that  all  of  the  above  analyses  assumes  that  the  fog  is 
d;lBtribvted  in  a nonodisperse  population,  tnd  we  have  initiated  work  to 
understated  the  imporbcuice  of  the  distribution  effects  on  these  measurements. 
For  instance,  we  oimuloaneously  stq>port  the  above  transmission  diaignostics 
with  scattering  measurements  (6328  %)  at  two  angles.  The  scattered  intensity 
ratio  for  two  angles  may  then  be  used  as  an  indicator  of  droplet  radius. 
However,  due  to  the  large  droplex  radius  to  wavelength  ratio  for  this 
dieignotic,  Mie  calculations  chow  very  complicated  scattering  patterns  develop- 
ing as  a function  of  radius  (Fig.  B-19).  Th.>  result  is  depicted  in  Fig. 

B-20  which  shows  a signi  fica’Ut  amount  of  ambiguity  for  the  scattered 
intensity  ratio  at  tvo  selected  angles.  Much  of  this  ambiguity  may  be 
removed  by  utilizing  measurements  at  angles  very  near  zero  degrees.  However, 
this  is  not  practical  for  ovir  experimental  arrangement.  Hence,  thiu  two 
angle  measurement  technique  is  not  a good  primary  indicator  of  droplet 
radius,  but  does  serve  as  a good  technique  for  confirmatior.  of  the  radius 
measured  via  the  transmission  diagnostics.  During  the  actual  experiments  we 
first  found  the  droplet  radius  using  the  transmission  data  and  checked  the 
Inferred  value  with  the  two  angle  scattering  predictions.  In  nearly  all 
cases  excellent  agreement  was  observed. 

The  in^jortance  of  this  fact  becomes  apparent  when  we  analyze  the  effect 
of  polydiapersion  or  this  backup  diagnostic  of  droplet  radius  measurements. 

To  do  this  we  have  developed  a computer  code  for  Mie  calculations  subject  to 
a polydlsperse  distribution  of  aerosols.  The  program  Incorporates  the  ZOLD 
distribution  re/iewed  in  Ref.  12.  The  v’idth  and  skewness  of  this  distribution 
are  specified  by  and  the  remaining  parameter  is  r^j,  the  modal  radius  of  the 
population.  For  nerrow  distributions  r * the  mean  radius  of  a norm^ 
distribution,  and  the  equiv8J.ent  width  of  a normal  distribution  is  o * 

Figvire  B-21  illustrates  the  density  fUi?ction  for  rjjj  ■ 2.5  pm  and  various 
Oq,  normalized  to  the  modal  value.  The  effect  of  this  distribution  with  a 
small  width  of  Oq  « 0.1  has  been  calculated  f^r  the  two  angle  scattering 
diagnostic.  Fig.  B-22.  This  figure  predicts  veiy  obvious  differences  in 
the  scattering  ratio  when  compared  with  the  monodisperse  results  of  Fig. 

B-20.  Furthermore,  our  mcasuremfents  ax’e  not  consistent  with  Fig.  B-22. 
Assioming  that  the  ZOLD  distribution  function  is  reasonable  for  modeling  cur 
problem,  the  measurements  would  infer  a real  distribution  significantly 
narrower  than  ■ 0.1.  However,  modeling  and  experiments  should  be  extended 
to  further  resolve  uncertainties  surrounding  the  distribution  dilemma 

To  accomplish  the  fog  characterization  measurements  discussed  above  we 
probe  about  2 m of  the  fog  with  coaxial  HeNe  eund  COp  lasers.  The  beams  were 
combined  and  directed  through  a ZnSe  window  into  the  propagation  cell.  After 
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passing  through  the  fog  aocl  exiting  through  another  window,  tht*  two  Lftair.i. 
are  aeparated  end  deflected  to  their  respective  detectors.  At  6328  8 
{chopped)  we  utilise  a simple  photodiode/operational  amplifier  system  to  pro- 
duce the  signal  for  phase  pensitive  detoctior  "by  a lock-iD-amplifier.  At 
10.6  microns  nufficient  sensitivity  for  signal  processing  was  obtained  from 
a Coherent  Radiation  Laboratories  Model  201  power  meter.  These  transmission 
signals  are  simultaneously  recorded  on  a strip  chart.  Concurrently  with  the 
transmission  recc-rdings,  we  record  scattered  light  signals  at  30  and  6o 
degrees  from  the  probing  HeNe  beam.  Then  during  data  analysis,  we  ccllect 
the  appropriate  signal  information  for  fog  characterization  recorded  just 
prior  to  the  time  when  the  high  power  beam  was  propagated  through  the  cell. 

A photograph  of  the  optical  system  is  provided  in  Fig.  B-23. 


5 . Spectrophone 

During  high  power  propagation  through  a fog,  some  of  the  intercepted 
energy  will  be  lost  to  scattering  and  evaporation  of  droplets,  and  some  will 
go  directly  into  heating  of  the  ambient  gas.  The  fraction  of  energy  t.hct  ie 
deposited  into  gas  heating  dictates  the  strength  of  the  thermal  blooming. 
Therefore,  it  is  important  to  measure  or  otherwise  document  this  amount  of 
heating.  One  very  sensiti  ve  method  of  achieving  this  information  .is  to 
measure  the  strength  of  pressure  waves  generated  from  the  gas  heating?"^ 

These  perturbations  can  be  detected  by  an  instrument  called  a spectrophone?5“30 
In  this  section  we  discuss  a spectrophone  developed  for  these  fog  propagation 
experiments. 

In  contrast  with  closed  cell  spectrophones  which  ust  Barocell^^ 
type  monoireters  as  sensors,  che  system  described  here  operates  with  a 
propagation  cell  that  is  not  sealed.  To  accon^lish  this,  the  laser  BOTu:*cft 
must  be  chopped.  The  spectrophone  electronics  are  then  "tuned”  to  the 
chopping  frequency  in  order  to  synchronously  detect  the  laser  generated 
pressure  waves.  Ibe  actual  detector  of  the  pressure  waves  is  a micro- 
phone? *30  jjj  system  that  we  have  constructed,  the  microphone 

is  followed  by  an  ultra-low  noise  preaa^lifier  (Fig.  B-24).  The  output  of 
this  preamplifier  is  filtered  by  a preselector  before  final  BWiasur iment  with 
a Icck-in-anqplifier. 


Table  B-VI  indicates  the  basis  characteristics  of  the  components  used  in 
our  spectrophone.  The  information  srjggesta  that  for  operating  conditions 
durir\g  which  no  background  noise  exists  the  minimum  detectable  pressure 
signal  will  be  dictated  by  the  noise  characteristics  of  the  microphone 
rather  than  the  subsequent  electronics.  In  practice,  however,  ambient 
laboratory  sound  levels  limit  the  minimum  detectable  pressure  signals.  This 
fact  ce.u  be  illustrated  by  the  follcifing  prelimiaary  apectrophone  tssts. 
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The  epectrophone  head  (consisting  of  the  microphone  and  the  praaiaplifler) 
vas  mounted  in  uhe  middle  of  a one  meter  propagation  cell.  The  cell  diame''’er 
was  about  3.6  cm.  Initial,  tests  were  conducted  with  germani  jr.  windowb  closing 
off  the  cell.  The  windows  were  used  only  to  hold  in  the  call*  ration  ges  and 
were  not  pressure-scaled.  A 6 watt  COp  beam  was  aimed  coaxially  down  the 
cell  and  chopped  at  500  Hz.  The  beam  radius  we.:  about  0.3  .:m.  Thr^  dimen- 
si  ojis  were  chosen  to  be  an  approximate  one-tentn  scaling  of  the  1 irger  fog 
propagation  experiment.  Tho  test  cell  wes  flooded  with  OOg,  and  the  Icck-in- 
amplifier  was  adjusted  for  maximum  response.  For  125  seconds  averaging  time 
and  an  effective  brnHwidth  of  0.001  Hz.;  the  output  vas  36  mV  nns.  Tl.e  beam 
wes  then  blocked  so  that  the  speetrophone  rospended  to  only  backgrouna  noise. 
The  background  signal  level  war  0.002  mV  rms.  These  measurements  imply  a 
signal  to  noise  ratio  of  1.?  x 10^.  Bince  Ihe  absorption  coefficient  of  CO2 
is  about  l.it  X 10“^  noise  equivalent  absorpti'fity  of  the  sptetro- 

pheas  is  about  7.6  x 10"*^  cm“^. 


When  the  cell  windows  vere  'removed  so  that  the  laboratory  background 
noise  had  its  full  impner  on  the  speetrophone j the  background  level  increased 
tc  0.035  itV  rms.  Hence,  the  signal  to  noise  ratio  for  the  completely  open 
cell  was  1.0  X 10^.  The  resultant  noise  equivalent  absorptivity  becomes 

1.'4  X 10"^ 


The  lQportdr.t  que.^tion  tc  be  answered  next,  is  whether  or  not  these 
spectrojjhone  characteristics  are  sufficient  for  the  full  scale  fog  propaga- 
tion exi.G.''ifflent.  References  25  Oiid  31  shew  |hat  the  laser  generated  pressure 
waves  will  have  to  be  proportional  to  uP/ita  where  a is  the  absorption 
coefficient  (cm”-^),  P is  the  laser  power  (W;  and  is  the  beam  radius  (cm). 
For  ova*  initia.l  speetrophone  tests,  this  q'.^antity  was  3.0  x lO”^.  In  the 
casfj  of  the  fog  propagation  experiment,  we  will  have  an  aerosol  absorption 
coefficient  of  about  1.0  x 10“^  cr.”^,  a power  of  5 kW  or  more  and  a beam 
radius  of  approximately  2.5  'm.  The  quantity  discussed  above  becomes  2.6  x 
10"^.  Tills  representE  nearly  an  order  of  magnitude  increase  in  the  signal 
to  noi»8  ratio  for  the  fog  experiments  and  pressures  still  within  the  dynamic 
range  of  the  microphone.  In  conclusion,  this  speetrophone  electronics  design 
is  satisfactory'  for  the  larger  fog  propagation  experiments. 

6.  Target  Plane  Diagnostics 


In  the  target  plane,  the  beam  la  attenuated  by  reflection  from  two  NaCl 
wedges  and  is  deflected  to  a (^canned,  linear  array  of  pyroelectric  detectors. 
Its  32  elem!int8  ccui  be  read  at  h 5 Wlz  rate,  but  only  500  Hz  was  utilized  in 
this  exper.iuent  - The  array  triggering  vas  synchronized  •*•0  the  chopper  in  the 
high  power  laser  oscillator  with  an  adjustable  time  delay  to  optimize  phasing. 
In  this  fashion,  ve  vert  able  to  neabui'e  a profile  of  the  beam  as  it  was 
translated  acrosa  the  ai’ruy.  Tee  array  '■•utput  wus  fed  into  an  adjustable 
gain  amplifier  and  then  to  an  operations,!  amplifier  follower  circuit  to 
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eliminate  loading  effects  on  the  detection  systems  calibration.  The  output 
signals  from  the  follower  were  recorded  on  a video  tap  recorder  used  as  a 
wide  bandwidth  PM  data  recorder.  These  signals  were  used  la^ur  to  reconstruct 
the  fog  induced  thermal  blooming  profiles. 

7.  Data  Analysis 

Analysis  of  the  recorded  beam  profiles  was  accomplished  by  digitizing 
the  array  generated  prollles  using  a Tektronix  Digital  Processing  Oscillo- 
scope. The  digitized  data  was  stored  on  floppy  discs  for  analysis  by  a PDP- 
11  minicomputer . Application  of  this  system  enabled  us  to  reconstruct  the  high 
power  beam's  two-dimensi  onal  intensity  distribution,  measure  measure 

the  increased  area  of  the  distorted  beam,  and  calculate  the  power  delivered 
to  the  target  by  integrating  over  the  intensity  distribution.  The  results 
of  this  analysis  are  presented  next. 

a.  RESULTS 

An  example  of  the  reconstructed  intensity  distributions  is  provided  in 
Pig.  B-25.  Here,  we  illustrate  on  an  arbitrary  intensity  scale  that  fog 
induced  thermal  blooming  (increased  beam  area)  is  reduced  merely  by  increasing 
the  laser  power  from  about  1.8  kW  to  Just  over  6.2  kW.  The  spatial  scale  is 
identical  for  both  runs.  The  area  of  the  lower  power  run  is  2 times  the 
lUidistorted  beam  area.  In  the  second  case,  the  beam  was  very  effective  in 
destroying  much  of  the  fog,  end  it  burned  its  way  to  the  target  with-  only  a 
2 percent  area  increase  over  the  xmdistorted  case.  This  reduced  thermal  dis- 
tortion was  achieved  even  though  oL  of  the  higher  power  run  was  U3  percent 
greater.  Referring  to  the  classical  model,  thermal  blooming  decreased  even 
t;hen  the  standard  nonlinear  distortion  parameter  increased  by  a factor  of 
about  5.  Ordinarily  for  molecular  absorption  this  would  have  resulted  in 
about  a 5-fold  increase  in  beam  area  going  from  the  lower  to  the  higher  power 
case,  but  instead  we  observed  a 2-fold  decrease  due  to  the  increased  burn- 
thrcugh  of  the  higher  power  case.  Our  attempts  to  comprehensively  analyze 


measurements  of  these  new  nonlinear  propagat.^on  effects  are  presented  in  the  I 

following  discussion  of  the  observed  parameters  ^rel*  ^rel'  I 
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Figure  B-26  sunmariaes  the  measuremeutB  of  for  the  range  of 

experiirental  conditions  listed  in  Table  1.  Irel  defined  as 

where  is  indicative  of  the  molecular  absorber  or  the  unmodified  fog 

conditions  in  the  cell.  The  sane  absorption  coefficient  is  used  deriving 
from  Eq.  (fl-T).  Figure  B-26  shows  that  our  experimental  system  measured 
molecular  absorption  (CO2  diluted  with  N2)  induced  thermal  distortion  that 
was  consistent  with  the  cleissical,  empirical  model  of  Eq. (B-6).  Furthermore, 
when  the  relative  hiaoidity  of  the  cell  was  forced  to  near  saturation  but  no 
fog  droplets  were  present,  we  again  measured  thermal  distortion  that  was  well 
described  by  the  molecular  absorption  model.  However,  in  this  case,  the 
dominate  absorption  coefficient  resulted  from  water  vapor. 

The  bresikdovn  of  the  empirical  model  is  dramatically  displayed  in  Fig. 
B-26  for  data  collected  when  fog  was  the  dominate  absorber  in  the  cell. 

During  these  conditions  P was  increased,  v was  decreased  or  was  increased 

to  force  larger  values  of  N^.  But  rather  than  observing  increased  thermal 
distortion,  we  measured  increasingly  better  propagation.  We  were  able  to 
deliver  up  to  10,000  times  more  peak  intensity  to  the  target  than  would  be 
expected  based  on  the  linear  absorption  due  to  the  fog  at  low  intensities. 

The  enhanced  propagation  resulted  from  forced  vaporization  of  the  droplets, 
enabling  Increased  transmission  and  decreased  thermal  distortion.  The 
following  discussion  will  enable  us  to  differentiate  between  the  two  effects. 

2-  Vel 

The  extent  of  decreased  thermal  distortion  is  illustated  in  Fig.  B-27 
which  shows  the  ratio  of  the  area  in  the  1/e  intensity  contour  of  the  undis- 
torted beam  to  the  distorted  profile.  The  data  is  plotted  as  a function  of 
the  same  used  in  Fig.  B-26.  The  empirical  model  is  used  in  this  case  to 
indicate  the  approximate  parametric  form  of  Apei  for  molecular  absorption 
induced  thermal  distortion.  The  data  shows  significantly  less  thermal 
blooming  than  expected.  However,  no  new  parametric  dependencies  on  Nf  sore 
observed  due  to  the  fact  that  many  different  P,  v,  and  fog  conditions  were 
utilized  dxirlng  these  measurements.  We  stress  at  this  point  that  based 
on  the  UTimodified  fog  characteristics  is  no  longer  a good  measure  of  para- 
metric effects. 


We  display  the  amount  of  enhanced  transmission  through  the  fog  in  Fig. 
B-23.  Pfei*  which  is  defined  as  the  power  traesmitted  to  the  target  plane 
through  the  1/e  intensity  contours  of  the  distorted  beam  in  proportion  to  the 
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undistorted  "beam,  is  given  as  a function  of  the  unmodified  optical  thickness 
of  the  fog,  ®g3ct^'  parameter  shoxild  illustrate  extinction  of  the  "beam 

independent  of  the  thermal  distortion  mechanism.  Our  measurements  show  that 
propagation  of  a beam  with  sufficient  power  to  force  droplet  vaporization 
provides  increased  transmission,  and  that  the  unmodifiei  fog  extinction 
coefficient  is  no  longer  useful  for  parameterization . 

4,  Spectrophone  Corrected 

During  some  of  our  data  runs,  we  successfully  operated  the  spectrophone 
in  conjunction  with  the  other  diagnostics.  The  spectrophone  measures  pressure 
waves  induced  in  the  ambient  gas  resulting  from  absorbed  energy  that  goes 
into  gas  heating  and  is  not  sensitive  to  that  absorbed  energy  which  is  used 
to  vaporize  particles.  Therefore,  the  effective  absorption  coefficient 
measured  by  the  spectrophone  ignores  the  unmodified  fog  parameters  and  responds 
directly  to  the  average  effective  absorption  coefficient  of  the  fog  as 
actively  modified  by  the  high  power  beam.  It  is  this  modified  fog  absorption' 
coefficient  which  causes  distortion.  The  spectrophone  measurements  are  cali> 
brated  by  comparing  signals  recorded  during  CO2,  molecular  absorption, 
blooming  experiments  with  signals  from  the  fog  experiments.  A simple, 
geometric  argument  is  used  to  account  for  blooming  induced  changes  in  beam 
radius  at  the  spectrophone  location.  The  effective,  modified  fog  absorption 
coefficient  is  utilized  to  calculate  new  values  of  and  Nf.  The  results 

are  compared  in  Fig.  B-29  with  the  same  parameters  based  on  the  unmodified  fog 
characterization.  The  spectrophone  corrected  measurements  agree  with  the 
empirical  model,  and  this  fact  supports  om*  code  calculations  that  show: 
thermal  blooming  induced  by  fog  agrees  with  the  empirical  model  if  the 
appropriate  modified  absorption  coefficient  within  the  beam  is  known. 

^ . Direct  Comparison  with  Code 

The  measured  and  values  were  compared  directly  with  the  modified 

nonlinear  propagation  code  taking  into  account  vaporization.  The  results 
are  shown  in  Figs.  B-30  and  B-31.  We  find  very  little  correlation  between 
the  measurements  and  the  code  predictions.  A similar  problem  was  encountered 
when  we  attempted  to  parameterize  the  measurements  by  the  (P/va)  modified 
absorption  coefficients  predicted  by  the  code  study.  The  scatter  is  too 
gi'eat  to  draw  any  useful  conclusions.  Without  further  extended  measurements 
and  code  exercises,  we  cannot  determine  whether  the  dominate  uncertainty  lies 
in  the  code  predictions  or  in  the  experimental  measurements. 


R77-922578-13 


SECTION  BIV 
CONCLUSIONS 


We  have  performed  a scaled  experimental  investigation  of  high  power 
laser  propagation  effects  in  fog.  The  beam  was  sufficiently  powerful  to 
vaporize  much  of  the  fog.  As  a result  we  observed  transmission  enhanced 
beyond  expected  levels  and  less  thermal  distortion  than  would  be  normally 
predicted.  The  standard  empirical  model  used  to  describe  thermal  blooming 
was  fovtnd  inadequate  when  based  on  the  standard  fog  absorption  coefficient. 

However,  if  a new  absorption  coefficient  was  identified  for  the  actively 
modified  fog  by  use  of  a spectrophone , then  the  empirical  model  aigreed  with 
the  measurements. 

In  support  of  the  measurements  we  have  developed  a nonlinear  propaga- 
tion code  that  takes  into  account  fog  vaporization  forced  by  the  high  power 
beam.  The  code  predicts  enhMced  transmission  and  decreased  thermal  distor- 
tion. Some  simplified  sealing  for  the  modified  fog  parameters  has  been 
identified,  and  as  in  the  experiments  when  the  modified  absorption  coefficient 
was  used  as  a parameter,  the  validity  of  the  empirical  model  was  retained. 

So  as  verified  by  spectrophone  measurements  and  detailed  code  calculations 
a usefxil,  simplified  parametric  scaling  law  has  been  identified  for  high 
power  propagation  through  fog.  To  describe  thermal  distortion  Eq.  B-33,  the 
beam  modified  fog  absorption  coefficient,  is  used  in  conjunction  with  Eqs.  B-6 
and  B^Tiwhile  to  predict  enhanced  transmission  Eq.  B-31  may  be  utilized  in  Eq.  B-32. 
Fineilly,  some  inconsistencies  were  found  for  a detailed  comparison  between 
the  absolute  code  predict ions and  the  experimental  measurements.  But  with 
the  studies  limited  to  this  report  we  have  not  been  able  to  identify  whether 
the  dominant  uncertainties  lie  in  the  new  code  or  in  the  experimental  mea- 
surements. Further  effort  should  be  extended  to  resolve  these  questions. 
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SECTION  BV 

SUGGESTIONS  FOR  FURTHER  WORK 


In  this  closing  section  we  list  items  that  should  he  considered  to  aid 
in  the  resolution  of  remaining  questions  with  respect  to  high  power  laser 
propagation  through  fog: 

1.  The  modified  code  should  he  exercised  over  a much  broader 
range  of  conditions.  This  would  aid  in  identification  of 
additional  simplified  scaling  possibilities,  and  identifi- 
cation of  real-scale  system  capabilities. 

2.  Additional  experiments  should  be  conducted  to  investigate 
limits  on  the  validity  of  the  modified  absorption  coefficient 
modeling  and  the  parametric  dependencies  of  aerosol  induced 
blooming . 

3.  Additional  fog  characterization  should  be  .nitiated  for  direct 
characterization  of  the  droplet  population  width  and  distribution. 

1*.  The  code  should  be  modified  to  account  for  fog  droplet 
populations  of  variable  width  and  distribution. 

5.  In  an  attempt  to  further  model  real-scale  meteorological 
conditions,  experiments  should  be  conducted  with  fogs  activated 
by  hygroscopic  nuclei. 

6.  Code  modifications  should  be  developed  to  predict  the  effects 
of  high  power  beams  on  fogs  with  hygroscopic  nucloi. 
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TABLE  I 


Contparison  of  experimental  parameters  and  approximate  field  conditions. 
P is  power;  a is  the  1/e  intensity  radius  of  the  bean;  L is  the  path- 
length;  I is  the  intensity,  v is  the  velocity,  aL  is  the  product  of 
the  absorption  coefficient  and  pathlength;  Nj.  Is  the  thermal  distortion 
scaling  parameter  for  a focused  beam  (defined  in  text);  r is  the  fog 
droplet  radius;  C is  the  fog  droplet* conren^'ati on;  and  m is  the  total 
water  mass  density  of  the  fog. 

.a’ . • 

i. 

LAB  FIELD 


Nondimensional 

Propagation 

Parameters 

F = 8.0 

aL  * 1-5 

Up  * 20-200 

Physical 

Propagation 

Parameters 


P/va 

P 


a 

L 


^2 

Ch2 

®j( Jitter) 


UO-200 
3 kW 
1.5  cm 
10  3 cm 
10-50  cm/s 
0 
0 


Fog  Parameters 


r 

C 

m 

visibility 

(6328  8) 


■ 2 X 10"**  cm 

■ 3-17  X 10**  car  3 

■ 1-5.7  gn/m3 

» 6.1*7  X 102  cm  - 3.21*  x 103 


8.3 

1-5 

20-200 


P/va  » 71-286 


2 X 105  cm 
50-200 
10“^^  n(-2/3 

10  prad 


2 X 10-**  cm 
1. 7-8.1*  X 102  cm-3 
0.0056-0.028  gm/m3 
1.3  X 105  cm  - 6.1*  X lo5 


Very  Dense  Generated  Fog 


Light  to  Medium 
Haturel  Fog 
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TABLE  II 

PANGE  OF  PARAMETERS  USED  IN  ACTUAL  EXPERIMENTS 

Power  = 0.89  - 6.2  kV 
Beam  radius  (l/e)  = 1.1*5  - 1.66 
Absorber  Length  = 1000  cm 
Length  to  detectors  » 1172  cm 
Beam  quality  * 1.76  - 1*.79 
Fresnel  number  = 2.6  - 7.9 
Velocity  =15-1+6  cm/s 

Fog  absorption  coefficient  (10.6  um)  * 0.71*  - 1+.26  x 10“^  cm~^ 


Fog  droplet  radius,  typical  * 2.0  ym 
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TABLE  III 


PROPAGATION  THROUGH  NONVOIaATILE  FOG 


Linear  Problem: 


Nonlinear  Problem; 


Clear  Atmosphere 


Linear  Problem: 


r = 2 X 10“^  cm 

Qext  * 3.07  X 10'^ 

C = 207  cm- 3 

“ext  (Eq.  1)  - 8.0  X 10“^  cnrl 
L ® 2 X lo5  cm 
» I/If,  (Eq.  2)  0.20 

r = 2 X 10"*'  cm 

Oats  = 2.5i*  X 10-1 
C * 207  cm- 3 

“abs  (Eq.  10)  » 6.6  x 10'^  cm-1 
F * 8.3 
V = 200  cm/s 
z = 2 X lo5  cm 
P/va  = 71.1*  w/cm2s 
Np  (Eq,  7)  » 25.7 

• Irel  * 0.0237 


“ext  » “abs  * 2 X 10-^  cnr^ 
L * 2 X 105  cm 

• I/Iq  (Eq.  2)  - 0.67 

“abs  * 2 ^ 10’^ 

F,  a,  V,  z P as  above 
Nj.  (Eq.  7)  « 10.1 

• ^rel 


Honlinep'  "’rob  1cm: 
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TABLE  IV 

PARAMETER  VABIATIOMS  OF  FIGS.  8-10 

Fixed:  = 8.75  X 

Z c 10  cm 
n ■ 1.0 

p » 1.18  X 10"^  cm‘^ 

c *1.0J/g°C 
P 

a = 1.^5  cm 
F - 1^.79 


Varied: 


Concentration 


k 

3.5  X 10 

U.O  X 


U.5  X lo"* 

5.9  X 10** 
k 

5.5  X 10 
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Table  v 

AEROSOL  GENERAJOR  DECISION  PROCESS 


^proximate  Geometric 
Standard  Deviation* 


of  Popxilation 

TVpe  of  Generator 

Comment 

< 1.05 

Berglxind  vibrating  orafice 

Too  few  per  second 

1.05  to  1.15 

Spinning  disc 

Too  few  per  second  and 
primary  droplets  too 
large 

•v  1,2 

Atomizer  with  impactor 
and  evaporative  control 

Droplets  too  small  and 
typically  used  with  oil 
rather  than  water 

1.2*  to  1.8 

Atomizer  with  in5)actor 

Droplets  too  small 

1.8  to  2.0 

Immersed  ultrasonic  transducer 

Too  few  per  second 

1.5  to  3.0 

Atomizer  with  acoustic  field 

Adequate  production 
rate,  size  range  and 
operating  pressures. 

« 

The  geometric  standard 
of  the  50  percent  point 

deviation  is  the  radius  of  the 
on  the  cumulative  distributions 

8U.1  point  to  the  radius 
of  the  droplet 

popiilatlon.  For  exa]i;>le,  if  the  geometric  etanderd  deviation  is  1.0  then 
all  droplets  have  the  same  size. 
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TABLE  VI 


SPECTROPHONE  COMPOMENTS 


Mlcr&phone 

Noise  (Total  A-wei^tel) 
(ji^ui  valent  1 NHz) 

Distortion  (3^  THD) 

Response  (He  IV/wB  for 
50*  2 ItHz) 

Dynamic  Range 


Electronic 


3C  db 


120  lo 


-70  dh  + 3 


90  db 


Absolute  Pres  pure 

fc.3  X 10“^  yB 

6.5  X 10”^  uB/./kz’ 


2.0  X 10  wB 


3.2  X IQ-^  V/wB 


3.2  X 10^ 


Preamplifier  HiM38lAJlj 

Eq,ui valent  Input  Noise 
(total,  10-10  ItHz) 
(at  1 kHz) 


5.5  nV/  ^ 
0.35  PA/  ^ 


1.6  X 10*^  PB 

5.1*  X 10*^  yB//Bz 

at  R,  * 1*7  kO 
L • 


DynaiAic  Range  (measured) 


Lock-In-/ 


Noise  (lypical  at  1 kHz) 


5 nV/»1iz 
0.03  pA/^Sz” 


1.6  z 10'^  pB/^Sz 

at  50  (1 


Bandvidth  (Selectable) 
Averaging  tine 
Accuracy 


0.001-100  Hz 
1.25  nS  - 125  S 

1% 


Chopper 

Frequency 
Bt ability 


a-  500  Hz 

0.1% 


Relative  to  2 x 10"  yB 
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FIG.  B-1 


76-10-143-3 
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FIC.  B-11 


I APPLICATION  OF  MODIFIED  SCALING  TO  CODE  PREDICTIONS 

EIVlPIRICAL  CAS  MODEL 

• CODE  PREDICTIONS  OF  r-,(J.  9 WITH  MODIFIED  SCALING 


LASER 


AEROSOL  PROPAGATION  CELL  SCHEMATIC 
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FIG.  B-16 


AEFIOSOL  PROPAGATION  CELL 


RL-77-154-D 
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DIAGNOSTICS  FOR  FOG  CHARACTERIZATION 


9 DATA 


EMPIRICAL  MODEL 


